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We derive a general effective action for quark matter at nonzero temperature and/or nonzero 
density. For this purpose, we distinguish irrelevant from relevant quark modes, as well as hard from 
soft gluon modes by introducing two separate cut-offs in momentum space, one for quarks, A 9 , and 
one for gluons, A s . We exactly integrate out irrelevant quark modes and hard gluon modes in the 
functional integral representation of the QCD partition function. Depending on the specific choice 
for A 9 and A 9 , the resulting effective action contains well-known effective actions for hot and/or 
dense quark matter, for instance the "Hard Thermal Loop" or the "Hard Dense Loop" action, as 
well as the high-density effective theory proposed by Hong and others. We then apply our effective 
^-jl ' action to review the calculation of the color-superconducting gap parameter to subleading order in 

. weak coupling, where the strong coupling constant g <C f . In this situation, relevant quark modes 

■ are those within a layer of thickness 2A 9 around the Fermi surface. The non-perturbative nature of 

£\j ' the gap equation invalidates naive attempts to estimate the importance of the various contributions 

via power counting on the level of the effective action. Nevertheless, once the gap equation has 
been derived within a particular many-body approximation scheme, the cut-offs A 9 , A 9 provide 
the means to rigorously power count different contributions to the gap equation. We recover the 
previous result for the QCD gap parameter for the choice A, ^ j/i < A, < n, where /i is the quark 
chemical potential. We also point out how to improve this result beyond subleading order in weak 
coupling. 
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PACS numbers: 12.38.Mh, 24.85.+p 



I. INTRODUCTION 



Quark matter at small temperature T and large quark chemical potential /i is a color superconductor Q, Q ■ While 
this discovery goes back to the late 1970's 3], wider interest in the phenomenon of color superconductivity has only 
recently been generated by the observation that, within a simple Nambu-Jona-Lasinio (NJL) - type model Q for 
the quark interaction, the color-superconducting gap parameter assumes values of the order of 100 MeV 0. Gap 
parameters of this magnitude would have important phenomenological consequences for the physics of compact stellar 
objects, and possibly even for heavy-ion collisions at laboratory energies of the order of ~ 10 AGeV. It is therefore 
of paramount importance to put the estimates from NJL-type models on solid ground and obtain a more reliable 
result for the magnitude of the gap parameter based on first principles. To this end, the color-superconducting gap 
parameter was also computed in quantum chromodynamics (QCD) \(i 1 71 l8L l9 UlCj . 

At zero temperature, T — 0, in weak coupling, g <C 1, and in the mean-field approximation, the gap equation for 
the color-superconducting gap parameter <p assumes the schematic form 



g 2 < 



C In 2 ( !=■ ) + (3 In 



(1) 



The solution is 



cf> = 2b(iexp^~j [I + 0(g)] . (2) 

The first term in Eq. is of leading order since, according to Eq. (J2J, g 2 In 2 (/i/ 4>) ~ 1. It originates from the 
exchange of almost static, long-range, Landau-damped magnetic gluons. One factor ln(/x/0) is the standard BCS 
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logarithm which arises when integrating over quasiparticle modes from the bottom to the surface of the Fermi sea, 
J dq/e q ~ In (///</>), where 



e q = ^{q - M )2 + <f (3) 

is the quasiparticle energy in a superconductor. The second factor ln(/i/</>) comes from a collinear enhancement 
~ ln(/Lt/e 9 ) in the exchange of almost static magnetic gluons. The coefficient £ determines the constant c in the 
exponent in Eq. J5J). As was first shown by Son 

C "7T (4) 

The second term in Eq. is of subleading order, g 2 ln(W0) ~ 5 -C 1. It originates from two sources. The first is the 
exchange of electric and non-static magnetic gluons 0,0 El El- In this case, the single factor ln(/i/0) is the standard 
BCS logarithm. The second source is the quark wave-function renormalization factor in dense quark matter |lllll2| . 
Here, the BCS logarithm does not arise, but the wave-function renormalization contains an additional ln(/i/e 9 ) which 
generates a ln(/x/</>). The coefficient j3 determines the prefactor b of the exponent in Eq. (J2J. For a two- flavor color 
superconductor, 

5/2 



256tt 4 (i— 2) exp( — ), (5) 



where Nf is the number of (massless) quark flavors participating in screening the gluon exchange interaction. The 
third term in Eq. is of sub- subleading order, ~ g 2 . The coefficient a determines the O(g) correction to the prefactor 
of the color-superconducting gap parameter in Eq. Since a has not yet been determined, the gap parameter can 
be reliably computed only in weak coupling, i.e., when the O(g) corrections to the prefactor are small. 

Due to asymptotic freedom the QCD coupling constant becomes small only at large momentum transfer. The 
typical momentum scale in dense quark matter is given by the quark Fermi momentum, kp = y/ i_i 2 — m 2 , where 
to is the quark mass. The Fermi momentum is equal to /1 up to terms of order 0(m 2 / '/i). Thus, j « 1 only for 
asymptotically large fi ^> Aqcd, where Aqcd is the QCD scale parameter. The range of \i values of phenomenological 
importance is, however, < 1 GeV. Although the quark density n is already quite large at such values of /i, n ~ 10 
times the nuclear matter ground state density, the coupling constant is still not very small, g ~ 1. It is therefore 
of interest to determine the coefficient of g in the O(g) corrections to the prefactor in Eq. J2J. If it turns out to be 
small, one gains more confidence in the extrapolation of the weak-coupling result J5J to chemical potentials of order 
~ 1 GeV. 

Let us mention that an extrapolation of the weak-coupling result @ for a two-flavor color superconductor, neglecting 
sub-subleading terms altogether and assuming the standard running of g with the chemical potential /i, yields values of 
4> of the order of ~ 10 MeV at chemical potentials of order ~ 1 GeV, cf. Ref. |2|. This is within one order of magnitude 
of the predictions based on NJL-type models and thus might lead one to conjecture that the true value of cf> will lie 
somewhere in the range ~ 10 — 100 MeV. However, in order to confirm this and to obtain a more reliable estimate of 
<f> at values of \i of relevance in nature, one ultimately has to compute all terms contributing to sub-subleading order. 

Although possible in principle, this task is prohibitively difficult within the standard solution of the QCD gap 
equation in weak coupling. So far, in the course of this solution terms contributing at leading and subleading order 
have been identified. However, up to date it remained unclear which terms one would have to keep at sub-subleading 
order. Moreover, additional contributions could in principle arise at any order from diagrams neglected in the mean- 
field approximation [ill Il3[ . Therefore, it would be ideal to have a computational scheme which allows one to 
determine a priori, i.e., at the outset of the calculation, which terms contribute to the gap equation at a given order. 

As a first step towards this goal, note that there are several scales in the problem. Besides the chemical potential 
there is the inverse gluon screening length which is of the order of the gluon mass parameter m g . At zero temperature 
and for Nf massless quark flavors, 

™*="/C£' (6) 

i.e., m g ~ g[i. Finally, there is the color-superconducting gap parameter 0, cf. Eq. (J2J). In weak coupling, g <C 1, 
these three scales are naturally ordered, <j> <C g/i -C /i- This ordering of scales implies that the modes near the Fermi 
surface, which participate in the formation of Cooper pairs and are therefore of primary relevance in the gap equation, 
can be considered to be independent of the detailed dynamics of the modes deep within the Fermi sea. This suggests 
that the most efficient way to compute properties such as the color-superconducting gap parameter is via an effective 
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theory for quark modes near the Fermi surface. Such an effective theory has been originally proposed by Hong |14l Il5| 
and was subsequently refined by others 0, 0, 0, 0] . 

At this point it is worth reviewing the standard approach to derive an effective theory [23, EI • In the most simple 
case, one has a single scalar field, <j), and a single momentum scale, A, which separates relevant modes, tp, from 
irrelevant modes, tp, <p = <p + ip. The relevant modes live on spatial scales L > 1/A, while the irrelevant modes live 
on scales 1 ^ 1/A « i. In the derivation of the effective action, one is supposed to integrate out the microscopic, 
irrelevant modes. Usually, however, this is not done explicitly. Instead, one constructs all possible operators Oi 
composed of powers of the field ip and its derivatives, which are consistent with the symmetries of the underlying 
theory, and writes the effective action as 

«W] = / $><0<(¥>). (7) 
Jx 4 

The coefficients, or vertices, gi determine the interactions of the relevant modes ip. A priori, they are unknown 
functions of the single scale A, g$ — <7i(A). All information about the microscopic scale I is contained in these vertices. 
Since the microscopic scale 1<L, the operators Oi are assumed to be local on the scale L. 

The effective action (JJJ contains infinitely many terms. In order to calculate physical observables within the effective 
theory, one has to truncate the expansion after a finite number of terms. One can determine the order of magnitude 
of various terms in the expansion Q via a dimensional scaling analysis which allows to classify the operators as 
relevant (they become increasingly more important as the scale L increases), marginal (they do not change under 
scale transformations), and irrelevant (they become increasingly less important as the scale L increases). To this 
end, one determines the naive scaling dimension of the fields, dhn(ip) = 6, from the free term in the effective action. 
Then, if the operator Oi consists of M fields ip and N derivatives, its scaling dimension is dim(Oi) = 6i = MS + N. 
The operator Oi is then of order ~ L~ Si . For dimensional reasons the constant coefficients gi must then be of order 
~ A d ~ Si , where d denotes the dimensionality of space-time. Including the integration over space-time, the terms in the 
expansion Q are then of order ~ (LA) d ~ Si . Consequently, relevant operators must have 6i < d, marginal operators 
Si = d, and irrelevant operators Si > d. At a given scale L, one has to take into account only relevant, or relevant 
and marginal, or all three types of operators, depending on the desired accuracy of the calculation. The final result 
still depends on A through the coefficients <7i(A). This dependence is eliminated by computing a physical observable 
in the effective theory and in the underlying microscopic theory, and matching the result at the scale A. 

There are, however, cases where this naive dimensional scaling analysis fails to identify the correct order of magni- 
tude, and thus the relevance, of terms contributing to the effective action. Let us mention three examples. For the first 
example, consider effective theories where, in contrast to the above assumption, the vertices gi are in fact non-local 
functions. Such theories are, for instance, given by the "Hard Thermal Loop" (HTL) or "Hard Dense Loop" (HDL) 
effective actions [22,[2j|. In these effective theories, valid at length scales L ~ l/(gT) or ~ l/(gn), respectively, there 
are terms g n A n in the effective action, which are constructed from a quark or gluon (or ghost) loop with n external 
gluon legs; A is the external gluon field with 6 = 1. The coefficients g n are non-local and do not only depend on the 
scale A < T, or < /i, but also on the relevant momentum scale 1/L ~ gT, or ~ g\i. Naively, one would expect g n 
to belong to a local n-gluon operator and to scale like A 4- ™. Instead, it scales like L n ~ 4 [23. For arbitrary n, the 
corresponding term g n A n in the effective action then scales like L A , independent of the number n of external gluon 
legs. 

The second example pertains to the situation when there is more than one single momentum scale A. As explained 
above, for a single scale A and a given length scale L, the naive dimensional scaling analysis unambiguously determines 
the order of magnitude of the terms in the expansion J7J). Now suppose that there are two scales, Ai and A2. Then, 
the vertices gi may no longer be functions of a single scale, say Ai, but could also depend on the ratio of A 2 /Ai. Two 
scenarios are possible: (a) two terms in the expansion J7J, say g n O n and g m O m , with the same scaling behavior may 
still be of a different order of magnitude, or (b) the two terms can have a different scaling behavior, but may still be 
of the same order of magnitude. In case (a), all that is required is that the operators O n and O m scale in the same 
manner, say L~ k , and that g n ~ A d ~ k , but g m ~ A d ~ k . If Ai -C A2, g m ^> g n , and thus the two terms are of different 
order of magnitude. In case (b), let us assume 1/L < Ai < A2, with A1/A2 ~ l/(Ai£) ~ e <C 1 and let us take the 
fields (p to have naive scaling dimension 6 = 1. Then, at a given length scale L, a term g„tp n , with a coefficient g n of 
order A d T n , can be of the same order of magnitude as a term g m ip m , m 7^ n, if the coefficient g m ~ A^ _m (A 2 /Ai) fc 
with k = d + m — 2n. Although the scaling behavior of the two terms is quite different as L increases, they can be of 
the same order of magnitude, if the interesting scale L happens to be ~ A2/A^. In both cases (a) and (b) the naive 
dimensional scaling analysis fails to correctly sort the operators Oi with respect to their order of magnitude. 

The third example where the naive dimensional scaling analysis fails concerns quantities which have to be calculated 
self-consistently. Such a quantity is, for instance, the color-superconducting gap parameter which is computed from a 
Dyson-Schwinger equation within a given many-body approximation scheme. In this case, the self-consistent solution 
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scheme leads to large logarithms, like the BCS logarithm in Eq. £[J. These logarithms cannot be identified a priori 
on the level of the effective action, but only emerge in the course of the calculation pj. 

In order to avoid these failures of the standard approach, in this paper we pursue a different venue to construct an 
effective theory. We introduce cut-offs in momentum space for quarks, A 9 , and gluons, A g . These cut-offs separate 
relevant from irrelevant quark modes and soft from hard gluon modes. We then explicitly integrate out irrelevant 
quark and hard gluon modes and derive a general effective action for hot and/or dense quark-gluon matter. One 
advantage of this approach is that we do not have to guess the form of the possible operators Oi consistent with 
the symmetries of the underlying theory. Instead, they are exactly derived from first principles. Simultaneously, the 
vertices gi are no longer unknown, but are completely determined. Moreover, in this way we construct all possible 
operators and thus do not run into the danger of missing a potentially important one. 

We shall show that the standard HTL and HDL effective actions are contained in our general effective action for a 
certain choice of the quark and gluon cut-offs A q , A g . Therefore, our approach naturally generates non-local terms in 
the effective action, including their correct scaling behavior which, as mentioned above, does not follow the rules of 
the naive dimensional scalin g an alysis. We also show that the action of the high-density effective theory derived by 
Hong and others [lj> [ljj Ua is a special case of our general effective action. In this case, relevant quark 

modes are located within a layer of width 2A q around the Fermi surface. 

The two cut-offs, A q and A 9 , introduced in our approach are in principle different, A q ^ A g . The situation is then as 
in the second example mentioned above, where the naive dimensional scaling analysis fails to unambiguously estimate 
the order of magnitude of the various terms in the effective action. Within the present approach, this problem does not 
occur, since all terms, which may occur in the effective action, are automatically generated and can be explicitly kept 
in the further consideration. We shall show that in order to produce the correct result for the color-superconducting 
gap parameter to subleading order in weak coupling, we have to demand A q < g\x <C A g < so that A q / A g ~ g <C 1. 
Only in this case, the dominant contribution to the QCD gap equation arises from almost static magnetic gluon 
exchange, while subleading contributions are due to electric and non-static magnetic gluon exchange. 

The color-superconducting gap parameter is computed from a Dyson-Schwinger equation for the quark propagator. 
In general, this equation corresponds to a self-consistent resummation of all one-particle irreducible (1PI) diagrams 
for the quark self-energy A particularly convenient way to derive Dyson-Schwinger equations is via the Cornwall- 
Jackiw-Tomboulis (CJT) formalism |24|. In this formalism, one constructs the set of all two-particle irreducible (2PI) 
vacuum diagrams from the vertices of a given tree-level action. The functional derivative of this set with respect to 
the full propagator then defines the 1PI self-energy entering the Dyson-Schwinger equation. Since it is technically not 
feasible to include all possible diagrams, and thus to solve the Dyson-Schwinger equation exactly, one has to resort to 
a many-body approximation scheme, which takes into account only particular classes of diagrams. The advantage of 
the CJT formalism is that such an approximation scheme is simply defined by a truncation of the set of 2PI diagrams. 
However, in principle there is no parameter which controls the accuracy of this truncation procedure. 

The standard QCD gap equation in mean-field approximation studied in Refs. [7|,|8|, |9| follows from this approach 
by including just the sunset-type diagram which is constructed from two quark-gluon vertices of the QCD tree-level 
action (see, for instance, Fig. 1181 below). We also employ the CJT formalism to derive the gap equation for the 
color-superconducting gap parameter. However, we construct all diagrams of sunset topology from the vertices of the 
general effective action derived in this work. The resulting gap equation is equivalent to the gap equation in QCD, 
and the result for the gap parameter to subleading order in weak coupling is identical to that in QCD, provided 
A 9 ^ gn <C A g < fi. The advantage of using the effective theory is that the appearance of the two scales A q and A g 
considerably facilitates the power counting of various contributions to the gap equation as compared to full QCD. 
We explicitly demonstrate this in the course of the calculation and suggest that, within this approach, it should be 
possible to identify the terms which contribute beyond subleading order to the gap equation. Of course, for a complete 
sub-subleading order result one cannot restrict oneself to the sunset diagram, but would have to investigate other 2PI 
diagrams as well. This again shows that an a priori estimate of the relevance of different contributions on the level 
of the effective action does not appear to be feasible for quantities which have to be computed self-consistently. 

This paper is organized as follows. In Sec. [n] we derive the general effective action by explicitly integrating out 
irrelevant quark and hard gluon modes. In Sec. IIIII we show that the well-known HTL/HDL effective action, as well 
as the high-density effective theory proposed by Hong and others, are special cases of this general effective action for 
particular choices of the quark and gluon cut-offs A g and A g , respectively. Section llVI contains the application of the 
general effective action to the computation of the color-superconducting gap parameter. In Sec. we conclude this 
work with a summary of the results and an outlook. 

Our units are Ti = c = ks = 1. 4- vectors are denoted by capital letters, — (fc ,k), with k being a 3- vector of 
modulus |k| = k and direction k = k/fc. For the summation over Lorentz indices, we use a notation familiar from 
Minkowski space, with metric g^ v — diag(+, — , — , — ), although we exclusively work in compact Euclidean space- 
time with volume V/T, where V is the 3- volume and T the temperature of the system. Space-time integrals are 

denoted as J^ T dr J y d 3 x = J x . Since space-time is compact, energy- momentum space is discretized, with sums 



5 



(T/V) J2 K = TJ2 n (l/V) J2 k . For a large 3-volume V, the sum over 3-momenta can be approximated by an integral, 
(1/V)^ k ~ J d 3 k/(27r) 3 . For bosons, the sum over n runs over the bosonic Matsubara frequencies lo^ = 2mrT, 
while for fermions, it runs over the fermionic Matsubara frequencies ui l n = (2n + l)irT. In our Minkowski- like 
notation for four-vectors, xq = t = —it, fco = —iuin^. The 4-dimensional delta-function is conveniently defined as 



II. DERIVING THE EFFECTIVE ACTION 



In this section, we derive a general effective action for hot and/or dense quark matter. We start from the QCD 
partition function in the functional integral representation fSec. lII A"|) . We first integrate out irrelevant fcrmion degrees 
of freedom fSec. HIB)) and then hard gluon degrees of freedom fSec. Ill Cj^l . The final result is Eq. (|53|l in Sec. Ill Dl We 
remark that the same result could have been obtained by first integrating out hard gluon modes, and then irrelevant 
fermion modes, but the intermediate steps leading to the final result are less transparent. 



A. Setting the stage 

The partition function for QCD in the absence of external sources reads 

Z = J VA exp{S A [A}} Z q [A] . 

Here the (gauge-fixed) gluon action is 



Sa[A] = f 
J x 



-\F^{X)F a (X) 



S gi [A] + S ghost [A} 



(8) 



(9) 



where = d^A® — d u A a ^ + g f abc A^A^ is the gluon field strength tensor, S g { is the gauge-fixing part, and Sghost the 
ghost part of the action. 

The partition function for quarks in the presence of gluon fields is 



Z q [A]= / XtyX>^exp{5,[4&V']} , 



where the quark action is 



S q [A, V] = / $[X) (iip x + M7o - m) tftX) 
Jx 



(10) 



(11) 



with the covariant derivative — — igA%(X)T a ; T a are the generators of the SU(N C ) C gauge group. In fermionic 
systems at nonzero density, it is advantageous to additionally introduce charge-conjugate fermionic degrees of freedom, 



Wl) = Cftl), i>c(X) = , ^(I)eC^(I) , $(X) = il%(X)C , (12) 

where C = i^jo is the charge-conjugation matrix, C _1 = = C T = —C, C~ 1 'y£C — — 7 M ; a superscript T denotes 
transposition. We may then rewrite the quark action in the form 

s q [A, *, *] = \ ^ g^(x, y) *(y) + lj x f^«(X) , (13) 

where we defined the Nambu-Gor'kov quark spinors 

9=( t) , 9 = $,$c) , (14) 



and the free inverse quark propagator in the Nambu-Gor'kov basis 

G -i (XY)= ([G+]-\X,Y) 



[G- Q ]-\X,Y)) ' (15) 
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with the free inverse propagator for quarks and charge- conjugate quarks 

{G±]- 1 (X,Y) = (ip x ±w -m)6W(X-Y) . (16) 
The quark-gluon vertex in the Nambu-Gor'kov basis is defined as 



ffi - ( 7 T° "t ) • (17) 



- 7 ^TJ 

As we shall derive the effective action in momentum space, we Fourier-transform all fields, as well as the free inverse 
quark propagator, 

= ^E e_iifX W' ( 18a ) 

§(X) = -j=£ e ^*(A), (18b) 
^(^F) - y^^e'^^^Q), (18c) 

= -±=Y,*- iPX A »a{P) ■ (18d) 

The normalization factors are chosen such that the Fourier-transformed fields are dimensionless quantities. The 
Fourier-transformed free inverse quark propagator is diagonal in momentum space, too, 

y i*>w - T y [Go]- 1 (A-) ) k 'Q ' [U) 



where [G^]" 1 ^) = # ± ^7o - m. 

Due to the relations (|12fl . the Fourier-transformed charge-conjugate quark fields are related to the original fields 
via ipc(K) = C l 4> T {— A), i/>c(A) = tp T (—K)C. The measure of the functional integration over quark fields can then 
be rewritten in the form 



>T>ip = JJ#(A)#(A) = Af J] di/j(K)dif>[K)dip{-K)dtp(-K) 



K (K-K) 

= Af' Yl d${K) d$(K) d$ c {K) d<ipc{K) = Af" f[ d^(K) dV(K) = X>§ , (20) 



(K, — K) (K-K) 



with the constant normalization factors A/", Af', Af". The last identity has to be considered as a definition for the 
expression on the right-hand side. 

Inserting Eqs. 1|18[) _ l|2(J|l into Eq. i|l(J|) . the partition function for quarks becomes 



Z g [A] = I V^V^f exp ~ § (g- 1 + ,g„4) * 



(21) 



Here, we employ a compact matrix notation, 

* {Go 1 + 9A) ^ = J2 [^'(^ Q) + 9A(K, Q)] #(Q) , (22) 



with the definition 



^(A,Q)EE- 7 i=f^(A-Q) . (23) 
The next step is to integrate out irrelevant quark modes. 
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B. Integrating out irrelevant quark modes 



Since we work in a finite volume V, the 3-momentum k is discretized. Let us for the moment also assume that there 
is an ultraviolet cut-off (such as in a lattice regularization) on the 3-momentum, i.e., the space of modes labelled by 
3-momentum has dimension D < oo. We define projection operators V\, V2 for relevant and irrelevant quark modes, 
respectively, 

^ 1 = T>i^> , ^2 = ^2* , *i = * 70P1 70 , *2 = * 70P270 ■ (24) 

The subspace of relevant quark modes has dimension N\ in the space of 3-momentum modes, the one for irrelevant 
modes dimension 7V 2 , with N\ + N 2 = D. 

At this point, it is instructive to give an explicit example for the projectors Vj .9. In the effective theory for cold, 
dense quark matter, which contains the high-density effective theory [Tj, [if! Ha, 03 LLa Oil discussed in Sec. IIII Bi as 
special case and which we shall apply in Sec. llVl to the computation of the gap parameter, the projectors are chosen 
as 

Vi(K, Q) = ( A * ^ ) 9(A 9 -\k- k F \) 4 4) Q , (25a) 

Vo(ko) = ( K + K®(\ k - k F\- A q ) \ s w (25h) 

'^ K ' Q > ~ { A+ + A k e(\k-k F \-A q ) ) dK >Q- (25b) 



Here, 



K ^2K [^ + e 7 o(7-k + m)] , (26) 



are projection operators onto states with positive (e = +) or negative (e = — ) energy, where E-^ = Vk 2 -I- m 2 is the 
relativistic single-particle energy. The momentum cut-off A q controls how many quark modes (with positive energy) are 
integrated out. Thus, all quark modes within a layer of width 2A q around the Fermi surface are considered as relevant, 
while all antiquark modes and quark modes outside this layer are considered as irrelevant. Note that, for the Nambu- 
Gor'kov components corresponding to charge-conjugate particles, the role of the projectors onto positive and negative 
energy states is reversed with respect to the Nambu-Gor'kov components corresponding to particles. The reason is 
that, loosely speaking, a particle is actually a charge-conjugate antiparticle. For a more rigorous proof compute, for 
instance, fe(Jf) = Cff(~K) using fa(-K) = ^(~K) 7o A+ k7o (for |fc - k F \ < A,) and 7o C[\±^ T C- X 7o = A k . 
In Sec. IIIII we shall discuss other choices for the projectors Vip, pertaining to other effective theories of hot and/or 
dense quark matter. The following discussion in this section, however, will be completely general and is not restricted 
to any particular choice for these projectors. 

Employing Eq. ipijl. the partition function l(2l1 becomes 

Z q [A]= f II ^»^«exp( I Y, M^ m ) ■ (27) 

•* n=l,2 V ra,m=l,2 / 

From now on, *&i,2) ^1,2 are considered as vectors restricted to the AT^-dimensional subspace of relevant /irrelevant 
3-momentum modes. The matrices Q^n> n = ■"•> 2) are defined as 

g~i(K, Q) - G^ nn {K, Q) + gA nn {K, Q) , (28) 

where the indices indicate that, for a given pair of quark energies ko, qo , the 3-momenta k, q belong to the subspace 
of relevant (n = 1) or irrelevant (n = 2) quark modes, i.e., Q~ x is an (N n x A^-dimensional matrix in 3-momentum 
space. The matrices Q^mi n ^ m, reduce to 

G^(K,Q)=gA nm (K,Q) , (29) 

since Q^ 1 is diagonal in 3-momentum space, i.e. Qo nm = for n ^ m. For a given pair of quark energies ko, qo, Qnm 
is a (N n x A r m )-dimensional matrix in 3-momentum space. 

The Grassmann integration over the irrelevant quark fields ^2, \&2 can be done exactly, if one redefines them such 
that the mixed terms ~ Qnmi n ^ m, are eliminated. To this end, substitute 



T = * 2 + G22 G21 *i , T = § 2 + #1 G12 G22 , 



(30) 
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where Q22 is the inverse of G 2 2 1 defined on the subspace of irrelevant quark modes. The result is 
Z q [A] = [ T&i 23*i exp \ ^ (Q^ - Q 22 Q^) t + i Tr 9 ln^ 1 



(31) 



The trace in the last term runs over all irrelevant quark momenta K, and not only over pairs (K, —K), as prescribed 
by the integration measure, Eq. 1121 III . This requires an additional factor 1/2 in front of the trace. A more intuitive 
way of saying this is that this factor accounts for the doubling of the quark degrees of freedom in the Nambu-Gor'kov 
basis. Of course, the trace runs not only over 4-momenta, but also over other quark indices, such as Nambu-Gor'kov, 
fundamental color, flavor, and Dirac indices. We indicated this by the subscript "g" . 
For a diagrammatic interpretation, it is advantageous to rewrite 



$11 - 012 022 Q 2 1 = 00,11 + 9& , 

where 

gB = gAn - gAi 2 G22 gA 2 i ■ 
The propagator for irrelevant quark modes, Q 22 , has an expansion in powers of g times the gluon field, 



G22 = 00,22 l)"g" [A22 00,22]" 

This expansion is graphically depicted in Fig. ^ 



(32) 
(33) 

(34) 



n=0 



+ 



+ 



+ 



FIG. 1: The full propagator for irrelevant quarks. The right-hand side symbolizes the expansion 13411 . The free irrelevant quark 
propagators So, 22 are denoted by double lines, the gluon fields A22 by curly lines. 



= 1 + j + j j I + j j'-'l 

FIG. 2: The diagrammatic symbol for the factor (1 + gAQa.22)^ 1 ■ 




FIG. 3: The term ^1 gB^i. A relevant quark field is denoted by a single solid line. 

Using this expansion, and suppressing the indices on A, Eq. Q33J1 can be symbolically written as 

gB = (1 + gAGo^r 1 gA , (35) 

which suggests the interpretation of the field Basa "modified" (non-local) gluon field. In the diagrams to be discussed 
below, the factor (1 + gAGo.22) will be denoted by the diagrammatical symbol shown in Fig. [3 With this symbol, 
the expression ^1 gB^i can be graphically depicted as shown in Fig. [21 
Since 

In g£ = In g J 2 - 9 n [00,22 ^22]" , (36) 

71=1 

the last term in the exponent in Eq. (|31|l also has a graphical interpretation, shown in Fig. 

This concludes the integration over irrelevant quark modes. Note that our treatment is (i) exact in the sense that 
no approximations have been made and (ii) completely general, since it is independent of the specific choice l|25|) for 
the projection operators. The next step is to integrate out hard gluon modes. 



o 




+ 






FIG. 4: The graphical representation of the term Tr g ln0 22 1 in Eq. 13111 . 



C. Integrating out hard gluon modes 

Combining Eqs. |J5J, (|31[1 . and the partition function of QCD for relevant quark modes and gluons reads 

3 = I V^ 1 V^ 1 VAexp{S[A,^ 1 ,^ 1 }} (37a) 

(37b) 



S[A, = S A [A] + \ *i + 5 S[A]} *i + \ Tr 9 ln0^[A] , 



where T>A = Y[ P dA(P). For the sake of clarity, we restored the functional dependence of the "modified" gluon field 
B and the inverse irrelevant quark propagator Q^ 2 on the gluon field A. 
The gluon action in momentum space is 

Sa[A] = -5 E Al{P l )[&-Xl{PuP2)Al{P 2 ) 

1 9 



E S pIp 2+ p 3 ,oK^(Pi,P2,Ps)A^P 1 )A^P 2 )A2(P 3 ) 



/•../•.../• 

2 



3! VVT3 

~3! (vlpj E ^U+ft+ft,oV^^(i\W(fli)^(ft)i45(P4) 
+Tr 9 , l lnW- 1 . (38) 
Here, A " 1 (Pi, P 2 ) is the gauge-fixed inverse free gluon propagator. To be specific, in general Coulomb gauge it reads 

[A r; b (a,*) - ^ [a„ r: wgu > ^ 

' (39b) 



where £c is the Coulomb gauge parameter and P M = (0,p). The vertex functions are 



V«%(P U P 2 , P 3 ) = - f abc [(P x - P 2 ) 7 .g Q(3 + (P 2 - P 3 ) Q 3/37 + (P 3 - PO0 5q7 ] , (40a) 

Vffi s = f abe r d {g ai g P s - gasgpy) + f ace f ebd {g aP g lS - g aS g^) + f ade f ebc (<m*M - g aJ 90s) ■ (40b) 



The last term in Eq. I|38|l is the trace of the logarithm of the Faddeev-Popov determinant, with the full inverse ghost 
propagator W _1 . The trace runs over ghost 4-momenta and adjoint color indices. 

Similar to the treatment of fermions in Sec. Ill Bl we now define projectors Qi, Q 2 for soft and hard gluon modes, 
respectively, 



where 



A X = Q X A , A 2 = Q 2 A, 



Q 1 (P 1 ,P 2 ) = 0(A g - Pl )5^ P . 
Q 2 (Pi,P 2 ) = e( Pl -A g )s^l P: 



(41) 

(42a) 
(42b) 



The gluon cut-off momentum A g defines which gluons are considered to be soft or hard, respectively. 
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We now insert A = A\ + A 2 into Eq. (|37() . The integration measure simply factorizes, VA = T)A\ T)A 2 . The action 
S[A, $1, ^1] can be sorted with respect to powers of the hard gluon field, 



= S[A 1 ,* 1 ,* 1 ]+A 2 J[A 1 ,* 1 ,* 1 ] --A 2 ^[Au^u^i] A 2 + Si[A 1 ,A 2 ,Vi,*i] 



(43) 



The first term in this expansion, containing no hard gluon fields at all, is simply the action (|37bjl , with A replaced by 
the relevant gluon field Ai. The second term, A 2 J , contains a single power of the hard gluon field, where 



J[A 1 ,V 1 ,y 1 
The first contribution, 



SA 2 



JB[-Al,*l,*l] + Jloopl^i] + JvlAi] 



A 2 =0 



2 V oA 2 J A Q 



(44) 



(45) 



arises from the coupling of the relevant fermions to the "modified" gluon field B, i.e., from the second term in Eq. 
(I37b|l . With the notation of Fig. [21 all diagrams corresponding to A 2 Jq can be summarized into a single one, cf. 
Fig. El It contains precisely two relevant fermion fields, ^1 and The second contribution, J\ OOVl arises from the 
terms Tr 9 lnC^ 1 and Tr g h lnW -1 in Eqs. I)37b(l . I)38|l. The loop consisting of irrelevant quark modes as internal lines, 
coupled to a single hard and arbitrarily many soft gluons, is shown in Fig. Finally, the third contribution, J7y, 
arises from the non-Abelian vertices, cf. Fig. [7| 



+ 



+ 



FIG. 5: The term A-^Jb- The hard gluon field is denoted by a dashed line, the soft gluon fields by wavy lines. 




FIG. 6: The fermionic contribution to the term AiJ\ oov . There is an additional contribution from ghosts with similar topology. 




+ 

FIG. 7: The term AiJv- 
The third term in Eq. (|43Jl is quadratic in A 2 , where 

5 2 SL4,*i,*i] 



A^[Ai,*i,*i] 



8A 2 SA 2 



n 22 [^i,*i,*i] 



(46) 



Here, A \ 2 is the free inverse propagator for hard gluons. Similar to the "current" J', cf. Eq. (|44|l . the "self-energy" 
II22 of hard gluons consists of three different contributions, 



n 22 [i4i,*i,*i] -n B [Ai,*i,*i] + nioop[-Ai] + n v [Ai] , 



(47) 
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, q + A+X 

FIG. 8: The term A2TI22A2 according to Eq. 14711 . The first diagram on the right-hand side corresponds to the term AhIIbAh. 
The second diagram is the fermion-loop contribution to A2TI100PA2; there is an analogous one from a ghost loop. The last two 
diagrams correspond to ^Ily^. 





FIG. 9: The term A 2 n B A 2 . 



which has a diagrammatic representation as shown in Fig. [5] The first two contributions on the right-hand side of 
Eq. 1(47(1 can be expanded as shown in Figs. I§l and 1 101 Figure ITT1 depicts the three- and four-gluon vertices contained 
in the last term in Eq. 1(47(1 . For further use, we explicitly give the first term, 



n B [Ai,*i,*i] 



*i 9 



5 2 B 
SA2SA2 



(48) 



A 2 =0 



Finally, we collect all terms with more than two hard gluon fields A2 in Eq. 1(43(1 in the "interaction action" for 
hard gluons, Si[A±, A2, $1, ^1]. We then perform the functional integration over the hard gluon fields ^2- Since 
functional integrals must be of Gaussian type in order to be exactly solvable, we resort to a method well-known from 
perturbation theory. We add the source term A2J2 to the action 1(37p(1 and may then replace the fields A 2 in Si by 
functional differentiation with respect to J2, at J2 = 0. We then move the factor exp{Si[Ai, S/6J2, ^1, ^i]} in front 
of the functional A 2 -integral. Then, this functional integral is Gaussian and can be readily performed (after a suitable 
shift of A2), with the result 



Z = / V^ 1 VHiVA 1 exp\ SLAi,*!,*!] 



1 



Tr g In A 22 



x exp < Si 



0J2 



exp 



(J + -h) A 22 {J + J2) 



(49) 



J 2 =0 



The trace over In runs over gluon 4-momenta, as well as adjoint color and Lorentz indices. We indicate this 
with a subscript "<j>" . Note that this result is still exact and completely general, since so far our manipulations of the 
partition function were independent of the specific choice IplSll for the projection operators Qi,2- The next step is to 
derive the tree-level action for the effective theory of relevant quark modes and soft gluons. 



D. Tree-level effective action 



In order to derive the tree-level effective action, we shall employ two approximations. The first is based on the 
principle assumption in the construction of any effective theory, namely that soft and hard modes are well separated 
in momentum space. Consequently, momentum conservation does not allow a hard gluon to couple to any (finite) 
number of soft gluons. Under this assumption, the diagrams generated by A2{J\ oop + J7v), cf. Fig. 00 will not occur 
in the effective theory. In the following, we shall therefore omit these terms, so that J = J$. Note that similar 
arguments cannot be applied to the diagrams generated by A 2 (ITi 00 p -fUy)^, cf. Fig.EB^H since now there are two 
hard gluon legs which take care of momentum conservation. 

Our second approximation is that in the "perturbative" expansion of the partition function 1)49(1 with respect to 
powers of the interaction action Si, we only take the first term, i.e., we approximate e Sl ~ 1. This is analogous to the 
derivation of the exact renormalization group in Ref. |25| , where it was shown that the corresponding diagrams are of 
higher order and can be neglected. In our case, diagrams generated by e Sl are those with more than one resummed 
hard gluon line. Even with the approximation e Sl ~ 1, Eq. 1(49(1 still contains diagrams with arbitrarily many bare 
hard gluon lines, arising from the expansion of 

In A^ 1 = In A - 22 - £ [ —^- (A , 22 n 22 )" , (50) 

n—1 
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FIG. 10: The fermionic contribution to the term A2Tl\ oop A2- 




FIG. 11: The term A 2 U V A 2 . 
and from the term Jb^-12Jb m Eq. (|49|l . when expanding 

oo 

A 22 = A , 22 £(-l) B ( n 22 A , 22 )" . (51) 

n=0 

With these approximations, the partition function reads 

Z = Jv^ 1 VWA X expjSeff^i, *!, , (52) 
where the effective action is defined as 

Settle, = SaIAJ + 1 *! + s B[Ai]} *i + ± Tr 9 ln^IAi] - ± Tr fl lnA^Ar, * lf ^] 

+ -J B [A U A 22 [Ai, Jh[A!, . (53) 

This is the desired action for the effective theory describing the interaction of relevant quark modes, "Fi, and soft 
gluons, A\. The functional dependence of the various terms on the right-hand side on the fields Ai, #1, \&i has been 
restored in order to facilitate the following discussion of all possible interaction vertices occurring in this effective 
theory. 




FIG. 12: The three- and four-gluon vertices in SUfAi], describing the self-interaction of soft gluons in Eq. I|53[l . 



The diagrams corresponding to these vertices are shown in Figs. 1121161 The three- and four-gluon vertices contained 
in iSa[Ai] are displayed in Fig. 1121 In addition, S^Ai] contains ghost loops with an arbitrary number of attached 
soft gluon legs. The topology is equivalent to that of the quark loops in Fig. 1 141 and is therefore not shown explicitly. 
The interaction between two relevant quarks and the "modified" soft gluon field, corresponding to <?2?[Ai] Wi, is 
depicted in Fig. El This is similar to Fig. except that now all gluon legs are soft. Diagrams where an arbitrary 
number of soft gluon legs is attached to an irrelevant quark loop are generated by Tr^lnC?^ 1 , cf. Fig. 1141 This is 
similar to Fig. 0] but now only soft gluon legs are attached to the fermion loop. The diagrams generated by the loop 
of a full hard gluon propagator, Tr 9 In A^ 1 , are shown in Fig. 1151 The first line in this figure features the generic 
expansion of this term according to Eq. (|5UII . where the hard gluon "self-energy" insertion II 22 , cf. Eq. (|47|) . is shown 
in Fig.|SJ The second line shows examples of diagrams generated by explicitly inserting II 22 in the generic expansion. 
Besides an arbitrary number of soft gluon legs, these diagrams also feature an arbitrary number of relevant quark 
legs. If there are only two relevant quark legs, but no soft gluon leg, one obtains the one-loop self-energy for relevant 
quarks, cf. the second diagram in the second line of Fig. 1151 The next two diagrams are obtained by adding a soft 
gluon leg, resulting in vertex corrections for the bare vertex between relevant quarks and soft gluons. The first of 
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FIG. 13: The term V 1 gB[A 1 ] $1 in the effective action CT . 



O^o + -0 + -0- + -Q:- 

FIG. 14: The term Tr q In [Ai] in the effective action 

these two diagrams arises from the n = 1 term in Eq. I|50(l . while the second originates from the n — 2 term. Four 
relevant quark legs and no soft gluon leg give rise to the scattering of two relevant quarks via exchange of two hard 
gluons, contained in the n = 2 term in Eq. (|5U|I . cf. the last diagram in Fig. El This diagram was also discussed in the 
context of the effective theory presented in Refs. [lH Il5| . cf. discussion in Sec. IIII 51 Finally, the "current-current" 
interaction mediated by a full hard gluon propagator, J$ A 2 2 Jb, Fig. El contains also a multitude of quark-gluon 
vertices. The simplest one is the first on the right-hand side in Fig. 1161 corresponding to scattering of two relevant 
fermions via exchange of a single hard gluon. 

The effective action l|53t is formally of the form J7J. The difference is that Eq. I|53[l contains more than one 
relevant field: besides relevant quarks there are also soft gluons. It is obvious that in this case there are many more 
possibilities to construct operators Oi which occur in the expansion (JJJ. As pointed out in the introduction, it is 
therefore advantageous to derive the effective action (|53[) by explicitly integrating out irrelevant quark and hard gluon 
modes, and not by simply guessing the form of the operators Oi, since then one is certain that one has constructed 
all possible operators occurring in the expansion (JJJ. 

As mentioned in the introduction, the standard approach to derive an effective theory, namely guessing the form 
of the operators Oi and performing a naive dimensional scaling analysis to estimate their order of magnitude, fails 
precisely when (a) there are non-local operators, or when (b) there is more than one momentum scale. Both (a) and 
(b) apply here. As we shall show below, the HTL/HDL effective action is one limiting case of Eq. I|53fl . and it is 
well known that this action is non-local. Moreover, as is obvious from the above derivation, there are indeed several 
momentum scales occurring in Eq. (|53|l . Let us focus on the case of zero temperature, T = 0, and, for the sake of 
simplicity, assume massless quarks, m = 0, fj, = kp. To be explicit, we employ the choice (|25|) for the projectors 
Vi.2 ■ In this case, the first momentum scale is defined by the Fermi energy /i. The propagator of antiquarks is 
~ l/(fco + /i + fc). If A q , A g < fi, the exchange of an antiquark can be approximated by a contact interaction with 
strength ~ l/fx, on the scale of the relevant quarks, L q ^> 1/A q > l//i, or of the soft gluons, L g >• 1/A a > 1/fx. 

The second momentum scale is defined by the quark cut-off momentum A q . The propagator of irrelevant quark 
modes is ~ l/(fco +fx — k). On the scale L q of the relevant quarks, not only the exchange of an antiquark, but also that 
of an irrelevant quark with momentum k satisfying \k — jj\ > A q is local, with strength ~ 1/A 9 . However, suppose 
that the quark cut-off scale happens to be much smaller than the chemical potential, A q -C fi. In this case, antiquark 
exchange is "much more localized" than the exchange of an irrelevant quark, l//i<C 1/A g . 

The third momentum scale is defined by the gluon cut-off momentum A g . The propagator of a hard gluon is 
~ 1/P 2 . On the scale L g of a soft gluon, the exchange of a hard gluon with momentum p > A g can be considered 
local, with strength ~ 1/A^. As we shall show below, in order to derive the value of the QCD gap parameter in weak 
coupling and to subleading order, the ordering of the scales turns out to be A q < g\x <C A g < fi. Thus, antiquark 
exchange happens on a length scale of the same order as hard gluon exchange, which in turn happens on a much 
smaller length scale than the exchange of an irrelevant quark, 1 / < 1/A S <C 1/A 9 . 

III. EXAMPLES OF EFFECTIVE THEORIES 

In this section we show that, for particular choices of the projectors Vi,2 in Eq. I|24|). several well-known, at first 
sight unrelated effective theories for hot and/or dense quark matter, are in fact nothing but special cases of the general 
effective theory defined by the action (|53|l . These are the HTL/HDL effective action for quarks and gluons, and the 
high-density effective theory for cold, dense quark matter. 
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FIG. 15: The term Tr 9 In AJ 2 [Ai, $1, $1] in the effective action 15311 . The first line corresponds to the generic expansion 1501 . 
with "self-energy" insertions IT22, as shown in Fig. [HI The second line contains some examples for diagrams generated when 
explicitly inserting the expression for II22- 



FIG. 16: The term Js^-2i,Ji3 in the effective action 15311 . The thick dashed line is a full hard ghion propagator, i.e., it has the 
expansion 1501 . The first diagram on the right-hand side of this figure results from the n = term of this expansion, while 
the next three diagrams originate from the n = 1 term. Even a single insertion of a hard gluon "self-energy" II22 gives rise 
to a variety of diagrams. Here, we only show the contributions corresponding to the first diagrams in Figs. El EH and the 
three-gluon vertex. The last diagram arises from the second term of the expansion shown in Fig. [5] 



A. HTL/HDL effective action 

Let us first focus on the HTL/HDL effective action. This action defines an effective theory for massless quarks 
and gluons with small momenta in a system at high temperature T (HTL), or large chemical potential /i (HDL). 
Consequently, the projectors Tip for quarks are given by 

Vi(K,Q) = Q(A q -k)5$ Q , (54a) 
V 2 (K,Q) = Q(k-A q )5$ Q , (54b) 

while the projectors for gluons are given by Eq. (|42|) . (We note that, strictly speaking, the quarks and gluons of 
the HTL/HDL effective action should also have small energies in real time. Since our effective action is defined in 
imaginary time, one should constrain the energy only at the end of a calculation, after analytically continuing the 
result to Minkowski space.) 

The essential assumption to derive the HTL/HDL effective action is that there is a single momentum scale, A q = 
A g = A, which separates hard modes with momenta ~ T, or ~ /i, from soft modes with momenta ~ gT, or ~ g/i. 
In the presence of an additional energy scale T, or /x, naive perturbation theory in terms of powers of the coupling 
constant fails. It was shown by Braaten and Pisarski [22( that, for the n-gluon scattering amplitude the one-loop term, 
where n soft gluon legs are attached to a quark or gluon loop, is as important as the tree-level diagram. The same 
holds for the scattering of n — 2 gluons and 2 quarks. At high T and small fi, the momenta of the quarks and gluons 
in the loop are of the order of the hard scale, ~ T. This gives rise to the name "Hard Thermal Loop" effective action, 
and allows to simplify the calculation of the respective diagrams. At large fi and small T, i.e., for the HDL effective 
action, the situation is somewhat more involved. As gluons do not have a Fermi surface, the only physical scale which 
determines the order of magnitude of a loop consisting exclusively of gluon propagators is the temperature. Therefore, 
at small T and large fj,, such pure gluon loops are negligible. On the other hand, the momenta of quarks in the loop 
are ~ /i. Thus, only loops with at least one quark line need to be considered in the HDL effective action. 
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In order to show that the HTL/HDL effective action is contained in the effective action l|53(l . we first note that 
a soft particle cannot become hard by interacting with another soft particle. This has the consequence that a soft 
quark cannot turn into a hard one by soft-gluon scattering. Therefore, 

gB[A 1 ] = gAu . (55) 

Another consequence is that the last term in Eq. Jb^-iiJb-, vanishes since J& is identical to a vertex between 

a soft quark and a hard gluon, which is kinematically forbidden. The resulting action then reads 

S larg0 7>L4i,*i,*i] = S A [A X \ + if i (g^ n + gAu) *i + i Tr q In g£ [Ax] - ~ Tr 3 In A^i, f 1; . (56) 

Using the expansion l|3(jf) we realize that the term Tr 9 XnQ^ generates all one-loop diagrams, where n soft gluon legs 
are attached to a hard quark loop. This is precisely the quark-loop contribution to the HTL/HDL effective action. 

For hard gluons with momentum ~ T or ~ fi, the free inverse gluon propagator is A^22 ~ T 2 or ~ /i 2 , while the 
contribution Hi oop to the hard gluon "self-energy" H47[) is at most of the order ~ g 2 T 2 or ~ g 2 ^ 2 - Consequently, ni oop 
can be neglected and H22 only contains tree- level diagrams, H22 = ng + Hy. Using the expansion l|50(l of Tr s In AJ 2 , 
the terms which contain only insertions of ny correspond to one-loop diagrams where n soft gluon legs are attached 
to a hard gluon loop. As was shown in Ref. p^. with the exception of the two-gluon amplitude, the loops with 
four-gluon vertices are suppressed. Neglecting these, we are precisely left with the pure gluon loop contribution to 
the HTL effective action. As discussed above, for the HDL effective action, this contribution is negligible. 

The "self-energy" n# contains only two soft quark legs attached to a hard quark propagator (via emission and 
absorption of hard gluons). Consequently, in the expansion (|50|l of Tr g InAj 2 , the terms which contain insertions of 
ny and Tig correspond to one-loop diagrams where an arbitrary number of soft quark and gluon legs is attached to 
the loop. It was shown in Ref. [22j that of these diagrams, only the ones with two soft quark legs and no four-gluon 
vertices are kinematically important and thus contribute to the HTL/HDL effective action. We have thus shown that 
this effective action, 5*htl/hdLi is contained in the effective action l|56|l . and constitutes its leading contribution, 

S'largcT/M = Shtl/hdl + hi g her orders . (57) 

For the sake of completeness, let us briefly comment on possible ghost contributions. Ghost loops arise from the term 
TrghlnW" 1 in SULAi]. Their topology and consequently their properties are completely analogous to those of the 
pure gluon loops discussed above. 

We conclude with a remark regarding the HDL effective action. According to Eq. 154(1 . at zero temperature and 
large chemical potential, a soft quark or antiquark has a momentum k ~ gfi, i.e., it lies at the bottom of the Fermi 
sea, or at the top of the Dirac sea, respectively. These modes are, however, not that important in degenerate Fermi 
systems, because it requires a large amount of energy fco ~ /1 to excite them. The truly relevant modes are quark 
modes with large momenta, k ~ fi, close to the Fermi surface, because it costs little energy to excite them. A 
physically reasonable effective theory for cold, dense quark matter should therefore feature no antiquark modes at all, 
and only quark modes near the Fermi surface. Such a theory will be discussed in the following. 



B. High-density effective theory 

An effective theory for high-density quark matter was first proposed by Hong 0] and was further refined by Schafer 
and others 0, 0, 0, 0] . In the construction of this effective theory, one first proceeds similar to our discussion in 
Sec.lrfland integrates out antiquark modes. (From a technical point of view, this is not done as in Sec.[H]by functional 
integration, but by employing the equations of motion for antiquarks. The result is equivalent.) On the other hand, 
at first all quark modes in the Fermi sea are considered as relevant. Consequently, in the notation of Sec. [H] the 
choice for the projectors Vi t 2 would be 

W,Q) = A ° k ) 4 4) q> (58a) 

W,Q) = ^ A °+) 4 4) q- (58b) 

Also, at first gluons arc not separated into soft and hard modes either. After this step, the partition function of the 
theory assumes the form ||SJ with Z q given by Eq. (|31|) . 

In the next step, one departs from the rigorous approach of integrating out modes, as done in Sec.[H] and follows 
the standard way of constructing an effective theory, as explained in the introduction. One focusses exclusively on 
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quark modes close to the Fermi surface as well as on soft gluons. However, since quark modes far from the Fermi 
surface and hard gluons are not explicitly integrated out, the effective action does not automatically contain the 
terms which reflect the influence of these modes on the relevant quark and soft gluon degrees of freedom. Instead, 
the corresponding terms have to be written down "by hand" and the effective vertices have to be determined via 
matching to the underlying microscopic theory, i.e., QCD. 

In order to further organize the terms occurring in the effective action, one covers the Fermi surface with "patches" . 
Each patch is labelled according to the local Fermi velocity, vjr = kfep/zx at its center. A patch is supposed to have a 
typical size A|| in radial (k) direction, and a size A_l tangential to the Fermi surface. The momentum of quark modes 
inside a patch is decomposed into a large component in the direction of vp, the particular Fermi velocity labelling 
the patch under consideration, and a small residual component, 1, residing exclusively inside the patch, 

k = (i\ F + l. (59) 

The residual component is further decomposed into a component pointing in radial direction, 1|| = vj?(vj? • 1), and 
the orthogonal one, tangential to the Fermi surface, lj_ = 1 — lii . The actual covering of the Fermi surface with such 
patches is not unique. One should, however, make sure that neighbouring patches do not overlap, in order to avoid 
double-counting of modes near the Fermi surface. In this case, the total number of patches on the Fermi surface is 

In the following, we shall show that the action of the high-density effective theory as discussed in Refs. [Til Il5l ITU ITU 
is contained in our effective action (|53|1 . To this end, however, we shall employ the choi ce 112511 and 14211 for the 
projectors for quark and gluon modes, and not Eq. 15%)) for the quark projectors. As in Refs. |l4ll5lll6Lll7llT8LIT^| . 
the quark mass will be set to zero, m = 0. We also have to clarify how the patches covering the Fermi surfaces 
introduced in Refs. [tH ITR flit ITtI ITsI IT^ | arise within our effective theory. It is obvious that the radial dimension Ai| 
of a patch is related to the quark cut-off A q . We simply choose A|| = A q . Similarly, since soft-gluon exchange is not 
supposed to move a fermion from a particular patch to another, the dimension A_l tangential to the Fermi surface 
must be related to the gluon cut-off A g . Again, we adhere to the most simple choice Aj_ = A g . Since A g < /j,, this 
is consistent with the matching procedure discussed in Ref. |l6j |. where the matching scale is chosen as Aj_ = 
(which is only slightly larger than The different scales A q , A g , and ijl are illustrated in Fig. El The modulus of 
the residual momentum 1 in Eq. I|59fl is constrained to I < max (A g , A g ). 

In Nambu-Gor'kov space, the leading, kinetic term in the Lagrangian of the high-density effective theory reads 

£ kin = ~ ^* 1 (A,vp) 7o ( iV Q D .y° Dc ) vMA,v F ) , (60) 

cf. for instance Eq. (1) of Ref. |l7|. Here, we have introduced the 4- vectors 

V = (l,v F ) , V" = (l,-v F ). (61) 

The covariant derivative for charge-conjugate fields is defined as = + igA^Tj . The contribution l|6U|l arises 
from the term ^ (GqIi + 9 An) in Eq. (|53|l . In order to see this, use V\ = V\ to write 



K,Q 



i^^^vphoi ( V Q L y° L ) *i(i,v F ). (62) 



In the last step, we have approximated k ~ /i + vp -1, which holds up to terms of order 0(l 2 /fj,), cf. Eq. I|59|l . This is a 
good approximation if the modulus of a typical residual quark momentum in the effective theory is I <ti max (A q , A g ) < 
fi. We have also introduced the 4- vector L M = (fco, 1) and, applying the decomposition l|59|) . we have written the sum 
over k as a double sum over vp and 1. The latter sum runs over all residual momenta 1 inside a given patch, while 
the former runs over all patches. With this decomposition, the spinors 'J'i are defined locally on a given patch 
(labelled by the Fermi velocity Vf), and depend on the 4- momentum L. Note that a Fourier transformation to 
coordinate space converts V ■ L — > iV ■ d. 

Now consider the term ^igA\\^\. Since An is not diagonal in momentum space, cf. Eq. (|23|l . in principle the 
two quark spinors \&i can belong to different patches. However, we have chosen the tangential dimension of a 
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FIG. 17: A particular patch covering the Fermi surface. The tangential dimension, Ax, is given by the maximum momentum 
transferred via a soft gluon, A H , while the radial dimension, Aii is defined by the maximum distance of relevant quark modes 
from the Fermi surface, A q . Also shown is a typical momentum transfer 1 via a soft gluon. 



patch such that a (typical) soft gluon can by definition never move a fermion across the border of a particular patch, 
|k — q| <C A g . Therefore, both spinors reside in the same patch and, to leading order, k ~ q ~ v F . With these 
assumptions we may write A~4 a (K -Q)A+ ~ V ■ A a {K - Q] ) 7o A£, A+4 a (K-Q)A^ ~ V ■ A a (K - Q) 7o A k . Then, 
introducing the residual momentum 1' corresponding to the quark 3-momentum q and defining L /AI = (go , 1') , the 
respective term in the effective action becomes 

vp,L,L' v 7 

In coordinate space, the sum of Eqs. (|62fl and l|63|) becomes Eq. ijfJUJl. 

Subleading terms of order 0(1/ fj,) in the high-density effective theory are of the form 

*o(x/,> - -l^iX^ho ± ( Dl " f 7 ! ^ _ Dlx _ i^j^j?* ) , (64) 

cf. Eq. (2) of Ref. [l7|. Here, = {0, (1 — vj?vj?) • D}, and similarly for D^ ± . The commutator of two gamma 
matrices is defined as usual, = (i/2)[y fJ ', 7"], and F^^ v T a = (i/g)[D± IJr , Dx. v \. As we shall see in the following, 
this contribution arises from the term — g 2 ^1 A12 Q22 A21 *i in Eq. (J53J. 

First, note that, with the projectors (|25(l . the irrelevant quark propagator C/22 contains quark as well as antiquark 
modes. In order to derive Eq. (j64(l . however, we have to discard the quark and keep only the antiquark modes. In 
essence, this is a consequence of the simpler choice for the projectors V\$, in the high-density effective theory of 
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Refs. [Tl 111 111 UtI HI Il9| . In this case, the propagator Q22 may be simplified. A calculation quite similar to that 
of Eqs. (|62|) and H63|) now leads (in coordinate space) to 

p—i 1 1 a (2n + iV-D \ , . 

Q22 =Qo,22 + . 9-422^70 r 3 (^ Q 2^-iV-Dc) ' (65) 

where T3 acts in Nambu-Gor'kov space. This result may be readily inverted to yield 

c If 1 f-iV-D \ n 



n=0 



Utilizing the projectors (J5SJ, one may also derive a simpler form for gA\2 and <?»42i- Consider, for instance, the term 
^1 gAu x &2- We follow the same steps that led to Eqs. I|62|). i.e., we assume that the spinors and ^2 reside in the 
same patch, such that k ~ q ~ vp. This allows to derive the identity jfi a (K — Q)A|F ~ K^J^]_{K — Q), where 
j^ml ^ ^ _ VfVf ) _ ^aj. Now introduce the 4-vectors , L /M , as in Eq. I|63[) . which leads to 

\^° A "**Ajk* S ^^(^V^ _ r °_ L , )raT )*2(L',v F ). (67) 

vf,L,L' 

We may add a term to the diagonal Nambu-Gor'kov components, which trivially vanishes between spinors "fi and 
\?2- This has the advantage that, in coordinate space, 

- ( 'V ilj! cx ) • ^ 

i.e., this term transforms covariantly under gauge transformations, and no longer as a gauge field. A similar calculation 
for 3^21 gives the result 5^21 = <7*4i2- Combining Eqs. ffify and JHEl, the term — g 2 ^1 A12 Q22 A21 ^1 corresponds 
to the following contribution in the Lagrangian, 

Taking only the n — term, and utilizing = g^ v — ia^", one arrives at Eq. 164|) . Note that our definition for 

transverse quantities, e.g. = {0, (I-vfVf)-A}, slightly differs from that of Refs. |14lll5l |. where A^_ = A^—V^V-A. 
However, both definitions agree when sandwiched between spinors ^i t 2 and ^2,1 • 

At order 0(l//x 2 ), besides the n = 1 term in Eq. I|69|) . there are also four-fermion interaction terms, cf. Eqs. (3-5) 
of Ref. H3- I n the effective action l|53|) . these contributions arise from the term Jb^-22Jb which originates from 
integrating out hard gluons. (Sin ce this is not done explicitly in the construction of the high-density effective theory 
in Refs. 0, 0, 0, 0> HE 0, this term is not automatically generated, but has to be added "by hand".) To 
leading order, this term corresponds to the exchange of a hard gluon between two quarks, cf. the first diagram on 
the right-hand side of Fig. ^j] If the quarks are close to the Fermi surface, the energy in the hard gluon propagator 
can be neglected, and Ao,22 ^ V-^-g- Since 1/A 2 , > 1 / , the contribution from hard-gluon exchange is of order 
0(1/ n 2 ). Four-fermion interactions also receive corrections at one-loop order, cf. Fig. 5 of Ref. |l5j. In Eq. (|53(l . they 
are contained in the term TrmA^~ 2 , see the last diagram in Fig. 1151 

Besides the quark terms in the Lagrangian of the high-density effective theory 0, 0, IrH 111 HI ] , there are 
also contributions from gluons. The first is the standard Yang-Mills Lagrangian — (1/4) F£ ll F£' / , cf. Eq. (1) of Ref. 
P3. This part is contained in the term in Eq. Ij53|) . cf. Eq. (J5J). The second contribution is a mass term for 

magnetic gluons, 

2 

771 

£ mg =-^fA a -A\ (70) 

cf. Eq. (19) of Ref. pfij . Eq. (18) of Ref. 0, or Eq. (27) of Ref. H3 > where m g is the gluon mass parameter {§J. 
This term has to be added "by hand" in order to obtain the correct value for the HDL gluon polarization tensor 
within the high-density effective theory. In Eq. (|53() this contribution arises from the n = 2 term of the expansion 1361) 
of TrlnC?^ 1 . The gluon polarization tensor has contributions from particle-hole and particle- antiparticle excitations. 
The latter give rise to C mg . While this term arises naturally within our derivation of the effective theory, it does 
not in the high-density effective theory of Refs. [li El 111 Il7l 111 HI] . because only antiquarks, but not irrelevant 
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quark modes, are explicitly integrated out. Irrelevant quark modes can then only be taken into account by adding 
the appropriate counter terms. 

Sometimes, the full HDL action is added to the Lagrangian of the high-density effective theory, cf. Eq. (8) of Ref. 
|17|. This procedure requires a word of caution. For instance, an important contribution to the HDL polarization 
tensor arises from particle-hole excitations around the Fermi surface. Such excitations are still relevant degrees of 
freedom in the effective theory. However, in order for them to appear in the gluon polarization tensor they would 
first have to be integrated out. Therefore, strictly speaking such contributions cannot occur in the tree-level effective 
action. Of course, in an effective theory one is free to add whatever contributions one deems necessary. However, one 
has to be careful to avoid double counting. As will be shown in Sec. II VI the full HDL polarization tensor will appear 
quite naturally in an approximate solution to the Schwinger-Dyson equation for the gluon propagator, however, not 
at tree-, but only at (one -)loop l evel. 

It was claimed in Refs. 0,13,0] that a consistent power-counting scheme within the high-density effective theory 
requires Aj_ = A||. In contrast, we shall show in Sec. IIVI that a computation of the gap parameter to subleading 
order requires A q = A\\ <C Aj_ = A g . This means that irrelevant quark modes become local on a scale l q 3> 1/A g , 
while antiquark modes become local already on a much smaller scale, lq ^> l//i, cf. discussion at the end of Sec. ITT1 
As mentioned in the introduction, for two different scales power counting of terms in the effective action becomes 
a non-trivial problem. While the high-density effective theory of Refs. 0, ITU H^ . UtI ITsl 0] contains effects from 
integrating out antiquarks, i.e., from the scale the effective action (|53H in addition keeps track of the influence 
of irrelevant quark modes, i.e., from physics on the scale 1/A 9 3> 1//U- Since all terms in the effective action (|53|l arc 
kept, one can be certain not to miss any important contribution just because the naive dimensional power-counting 
scheme is invalidated by the occurrence of two vastly different length scales. 

IV. CALCULATION OF THE QCD GAP PARAMETER 

In this section, we demonstrate how the effective theory derived in Sec. |H] can be applied to compute the gap 
parameter of color-superconducting quark matter to subleading order. For the sake of definitcncss, we shall consider 
a spin-zero, two-flavor color superconductor. 

A. CJT formalism for the effective theory 

The gap parameter in superconducting systems is not accessible by means of perturbation theory; one has to apply 
non-perturbative, self-consistent, many-body resummation techniques to calculate it. For this purpose, it is convenient 
to employ the CJT formalism j2J]. The first step is to add source terms to the effective action l|53l) . 

S eS [A x , *!, *i] — » 5 off [Ai, #i, tf i] + J X A X + X - A x KxAx + ^ (tf + Jf x tf i + , (71) 

where we employed the compact matrix notation defined in Eq. (|22|l . J\, Hi, and Hx are local source terms for 
the soft gluon and relevant quark fields, respectively, while K\ and /Ci are bilocal source terms. The bilocal source 
/Ci for quarks is also a matrix in Nambu-Gor'kov space. Its diagonal components are source terms which couple 
quarks to antiquarks, while its off-diagonal components couple quarks to quarks. The latter have to be introduced 
for systems which can become superconducting, i.e., where the ground state has a non- vanishing diquark expectation 
value, (ipxipi) 7^ 0. 

One then performs a Legendre transformation with respect to all sources and arrives at the CJT effective action 
[HI23 

r[A,%%A,g] = S c ff [A, — - Tr g In A -1 - i Tr g (D _1 A - l) 

+ ^Tr q hig- 1 + ^Tr 9 (G- 1 g-l)+T2[A,9,9,A,g] . (72) 

Here, S c s[A, ^f] is the tree-level action defined in Eq. (|53JI . which now depends on the expectation values A = (Ax), 
^ = (^1), and if? = (if>x) f° r the one-point functions of soft gluon and relevant quark fields. In a slight abuse of 
notation, we use the same symbols for the expectation values as for the original fields, prior to integrating out modes. 
This should not lead to confusion, as the original fields no longer occur in any of the following expressions. 

The quantities D^ 1 and G _1 in Eq. (JEJ are the inverse tree-level propagators for soft gluons and relevant quarks, 
respectively, which are determined from the effective action S c g, see below. The quantities A and Q are the expectation 
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values for the two-point functions, i.e., the full propagators, of soft gluons and relevant quarks. The functional 
is the sum of all two-particle irreducible (2PI) diagrams. These diagrams are vacuum diagrams, i.e., they have no 
external legs. They are constructed from the vertices defined by the interaction part of S e g, linked by full propagators 
A, Q. The expectation values for the one- and two-point functions of the theory are determined from the stationarity 
conditions 

The first condition yields the Yang-Mills equation for the expectation value A of the soft gluon field. The second 
and third condition correspond to the Dirac equation for and respectively. The effective action (|53|) contains 
a multitude of terms which depend on A, "J, and thus the Yang-Mills and Dirac equations are rather complex, 
wherefore we refrain from explicitly presenting them here. Nevertheless, for the Dirac equation the solution is trivial, 
since $1, <J/i are Grassmann-valued fields, and their expectation values must vanish identically, ^ = (^>i) = \f r = 
(^>i) = 0. On the other hand, for the Yang-Mills equation, the solution A is in general non-zero but, at least for 
the two-flavor color superconductor considered here, it was shown |2^,|2^| to be parametrically small, A ~ cf> 2 / (g 2 fj,) , 
where 4> is the color-superconducting gap parameter. Therefore, to subleading order in the gap equation it can be 
neglected. 

The fourth and fifth condition i|73|) are Dyson-Schwinger equations for the soft gluon and relevant quark propagator, 
respectively, 

A" 1 = D- 1 +11, (74a) 
g- 1 = G- l +V, (74b) 



where 



H = -2 ^ , (75a) 



E ee 2^ (75b) 



are the gluon and quark self-energies, respectively. The Dyson-Schwinger equation for the relevant quark propagator 
is a 2 x 2 matrix equation in Nambu-Gor'kov space, 

^-( K TVi-») + (S£)- (76) 

where E + is the regular self-energy for quarks and E~ the corresponding one for charge-conjugate quarks. The off- 
diagonal self-energies < & ± , the so-called gap matrices, connect regular with charge- conjugate quark degrees of freedom. 
A non-zero < & ± corresponds to the condensation of quark Cooper pairs. Only two of the four components of this 
matrix equation are independent, say [G + ] _1 + E + and <£>+, the other two can be obtained via [G~] _1 + E~ = 
CUG+}- 1 + E+} T C-\ $" ee 7o[$ + ] t 7o. Equation JUJ) can be formally solved for Q [13, 



g+ h- 
s+ g- 



(77) 



where 



g ± EE{[G ± ]- 1 + E ± -$T([G T ]- 1 + E T ) 1 $ ± } 1 (78) 

is the propagator describing normal propagation of quasiparticles and their charge- conjugate counterpart, while 

~± = -({G^}- 1 + ^y 1 $±g± (79) 

describes anomalous propagation of quasiparticles, which is possible if the ground state is a color-superconducting 
quark-quark condensate, for details, see Ref. 0. 
The tree-level gluon propagator is defined as 

D = 6ASA ■ (80) 
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Since we ultimately evaluate the tree-level propagator at the stationary point of T, Eq. (|73() . where ^ = ^ = 0, we 
may omit all terms in S e g, Eq. l|53ll . which are proportional to the quark fields. The only terms which contribute to 
the tree-level gluon propagator are therefore 

D ~* = ~JZ5A { Sa + \ Tlq ln a » " \ TTg ln A » ) ' (81) 
Using the expansions 1)34(1 . 1(36(1 . I|5()|l. and l|51|l. and exploiting the cyclic property of the trace, one finds 

= -^-2 Tr ^^^J + 2 Tl H^F7T + A22 MMj ■ (82) 

In order to proceed, note that the Dyson-Schwinger equations 1(7411 are evaluated at the stationary point of the effective 
action, where $ = W = 0, A ~ 0. For A = 0, the first term yields the free inverse propagator for soft gluons, Aq~ 1;l , 
cf. Eq. (|38|) . plus a contribution from the Faddeev-Popov determinant, (S 2 TT g } l lnW -1 /<5 ASA)a=o- The contributions 
from the three- and four-gluon vertex vanish for A = 0. Furthermore, according to Eq. ((23|1 . 

M 22 (A',Q) _ 1 - (4) =Y(KO-P) (83) 

This is a matrix in fundamental color, flavor, and Nambu-Gor'kov space, as well as in the space of quark 4-momenta 
K,Q. It is a vector in Minkowski and adjoint color space (F carries a Lorentz- vector and a gluon color index), as 
well as in the space of gluon 4-momenta P. We evaluate (5G 22 /SA)a=q using the expansion 1(34(1 . Only the term 
for n — 1 survives when taking A = 0. For ^ = $ = 0, we have Hb — 0, cf. Fig. and we only need to consider 
II 22 = rii 00 p + Hy. Then, the term = (SHy /8A)a=o corresponds to a triple-gluon vertex, cf. Fig. EI where two 
hard gluons couple to one soft gluon. The term (5Hi oop / 5 A) a=o is a correction to this vertex: it couples two hard gluons 
to a soft one through an (irrelevant) quark loop, cf. Fig. 1101 According to arguments well-known from the HTL/HDL 
effective theory, this vertex correction can never be of the same order as the tree-level vertex V*- 3 ^, since the two 
incoming gluons are hard. We therefore neglect (5Hi oc ,p/5A)a=o in the following. Similarly, V' 4 ' = (5 2 Uv/6A5A)a=o 
is a four-gluon vertex, cf. Fig. ^2 where two hard gluons couple to two soft ones, and (5 2 H\ oop /SA6A)a=o is the 
one- (quark-) loop correction to this vertex, cf. Fig. 1101 Applying the same arguments as above, we only keep V^- 4 \ 
Arguments from the HTL/HDL effective theory also tell us that to leading order we may approximate A 22 ~ Ao. 22 . 
Finally, utilizing the same arguments we approximate <5A 22 /<5A ~ — Ao j22 V^ Ao j22 . Then, the inverse tree-level gluon 
propagator of Eq. (|82() becomes 

D- 1 = A^ + ^ Tr 9 (g Q , 22 f G , 22 f ) - \ Tr s (a , 22 A , 22 V( 3 >) + \ Tr g (a , 22 V^) - 

U= (84) 

The second term represents an (irrelevant) quark-loop, while the third term is a hard gluon loop. The fourth term is a 
hard gluon tadpole. Finally, the last term in Eq. H84JI corresponds to a ghost loop necessary to cancel loop contributions 
from unphysical gluon degrees of freedom. Note that, in the effective theory, loop contributions involving irrelevant 
quarks and hard gluons occur already in the tree- level action ((53(1 . Therefore, such loops also arise in the inverse tree- 
level propagator 1)84(1 for the soft gluons of the effective theory. For the projection operators 1(42(1 and 154J1 the inverse 
tree-level propagator (|84|l is precisely the HTL/HDL-rcsummcd inverse gluon propagator. For small temperatures, 
T <C ;it, the contribution from the gluon and ghost loops is negligible as compared to that from the quark loop, 

2 

D- 1 - A^ + ^-H:, (00,22 f 0o,22 f) • (85) 
The inverse tree-level quark propagator is defined as 

G-^-2 fi8 «[A,*,*] m (86) 

For ^ = = 0, the last term in Eq. I(53|l does not contribute to G _1 , because it has at least four external quark legs, 
and the two functional derivatives S/Sfy amputate only two of them. The first and the third term in Eq. 1(53(1 

do not depend on ^i,^ at all, therefore 

A 2 Tr In A -1 

^=G&+9B[A] + i*$^ ■ (87) 
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Using the expansion formulae (|50|l and (|51() and the fact that H22 depends on "J, W only through rig, we obtain 

S^Sty y \ Sty Sty J 

We have exploited the fact that this expression is evaluated at ty = ty = 0, i.e., terms with external quark legs will 
eventually vanish. The trace runs only over adjoint colors, Lorentz indices, and (hard) gluon 4-momenta. Since A22 
is a hard gluon propagator, the contribution from IT22 to A22 may be neglected to the order we are computing, and 
we may set A 22 ^ A 0j2 2- Furthermore, (S 2 U B /StySty) A =o = -g 2 ^ Qo,22^, cf. Fig. El At ty = * = 0, A ~ we are 
left with 

G~ l = Go,\x ~ 9 2 Tr s (A 0)22 f £0,22 f) . (89) 

As was the case for the tree-level gluon propagator, also the tree-level quark propagator receives a loop contribution; 
here it arises from a loop involving an irrelevant quark and a hard gluon line. The term TQq^F under the gluon trace 
remains a matrix in the quark indices, i.e., fundamental color, flavor, Dirac, and quark 4-momenta. 

We now proceed to solve the Dyson-Schwinger equations l|74|l for the soft gluon and relevant quark propagator. To 
this end, we have to determine T 2 . Of course, it is not feasible to consider all possible 2PI diagrams. The advantage 
of the CJT formalism is that any truncation of T2 defines a meaningful, self-consistent many-body approximation for 
which one can solve the Dyson-Schwinger equations i|74|) . ln our truncation of T2 we only take into account 2-loop 
diagrams which are 2PI with respect to the soft gluon and relevant quark propagators A, Q, 

^2 = ~^Tr q , g (gfgf A) - y Tr,, 9 (^f £ ,2 2 f A) _^Tr 9 , g (gf£f A , 22 ) . (90) 

The traces now run over quark as well as over gluon indices. Consider, for instance, the term Q T Q T. It is a matrix 
in the space of fundamental color, flavor, Dirac and quark 4-momenta, of which the trace is taken through Tr g . In 
addition, due to the two factors T it carries two Lorentz-vector, adjoint-color, and gluon-4-momenta indices. The 
trace Tr g contracts these indices with the corresponding ones from the gluon propagator A. 

r 2 = Q + 2 Q, + Q 

FIG. 18: Diagrammatic representation of V2, Eq. I9LIH . 



The diagrams corresponding to Eq. (|90|l are shown in Fig. ^] The first two terms are constructed from the quark- 
gluon coupling ~ ty gB \f r . Using Eq. , one may either obtain an ordinary quark-gluon vertex ~ g ty A ty , involving 
one soft gluon and two relevant quark legs, or a vertex ~ g 2 ty AG22 Aty , with (at least) two soft gluon legs and two 
relevant quark legs. To lowest order, we approximate Q22 — G0.22, which neglects vertices with more than two soft 
gluon legs. Taking two ordinary quark-gluon vertices and tying them together to obtain a 2PI 2-loop diagram, we 
arrive at the first term in Eq. (|90|l . or the first diagram in Fig. ^| Taking one of the two-gluon-two-quark vertices 
and tying the legs together, one obtains the second term in Eq. H9U|) . or the second diagram in Fig. respectively. 
Finally, the third term/diagram arises from the last term in Eq. I|53l) . To lowest order, this corresponds to a four-quark 
vertex ~ g 2 ty T ty Ao.22^ T ty. Tying the quark legs together to form a 2PI diagram, one obtains the corresponding 
term/diagram in Eq. P0)l /Fig. ITSl 

The combinatorial factors in front of the various terms in Eq. H9U|) are explained as follows. In the first diagram, 
there are two ordinary quark-gluon vertices. According to Eq. I)53|l . each comes with a factor 1/2. Moreover, since 
there are two vertices, the diagram is, in the perturbative sense, a diagram of second order, which causes an additional 
factor 1/2 [31]. Finally, there are two possibilities to connect the quark lines between the two vertices. In total, we 
then have a prefactor — (1/2) 2 x 1/2 x 2 = — 1/4, where the minus sign arises from the fermion loop. The second 
diagram arises from the two-quark- two-gluon vertex, which already comes with a prefactor —1/2 in Eq. (|53|l . It is 
perturbatively of first order, and there is only one possibility to tie the quark and gluon lines together, so there is no 
additional combinatorial factor (and no additional minus sign) for this diagram. Finally, the third diagram arises from 
the four-quark vertex, (1/2) 1 7bAo, 22i /b, in Eq. (|53|l . This vertex comes with a factor 1/2 and is perturbatively of first 
order. However, there are two additional factors 1/2 residing in since Jb ~ (l/2)^ , r^> ', cf. Eq. Again, there 

are two possibilities to tie the quark lines together, so that, in total, we have a prefactor -1/2 x (1/2) 2 x2 = -1/4, 
where the minus sign again stands for the quark loop. 



23 




FIG. 19: Diagrammatic representation of Tj , after decomposing quark lines into relevant and irrelevant, as well as gluon 
propagators into soft and hard contributions. 

At this point, it is instructive to compare T 2 , Eq. (HP, in the effective theory with r^ CD which one would have 
written down in QCD at the same loop level. If D would be equivalent to the first diagram of Fig. [^J but now the 
quark and gluon lines represent the full propagators for all momentum modes, relevant and irrelevant as well as soft 
and hard. In order to compare with T 2 of the effective theory, we decompose the quark propagators into relevant 
and irrelevant modes, and the gluon propagator into soft and hard modes. One obtains the six diagrams shown in 
Fig. 1191 The first three are precisely the same that occur in r 2 of the effective theory, including the combinatorial 
prcfactors. The last three diagrams do not occur in T 2 of the effective theory, because they are not 2PI with respect 
to the relevant quark propagator Q and the soft gluon propagator A. Nevertheless, they are still included in the 
CJT effective action of the effective theory, Eq. (I72f) : opening the relevant quark line of the fourth diagram, wc 
recognize the loop contribution to the tree- level quark propagator G _1 , cf. Eq. i|89|) . Now consider the fifth term in 
Eq. (|72(l : here, this loop contribution to G _1 is multiplied with Q and traced over, which yields the fourth diagram in 
r^ CD . Similarly, opening the soft gluon line of the fifth diagram, we identify this diagram as the irrelevant quark- loop 
contribution to the tree-level gluon propagator -D -1 , cf. Eq. (|85(l . The third term in Eq. (|72J) . where this contribution 
is multiplied by A and traced over, then yields the fifth diagram of T^ D . Finally, the sixth diagram resides in the 
term ~ Tr g In A^ 1 of the tree-level effective action S e g, cf. Fig.^| Therefore, in principle, the CJT effective action 
(|72|l for the effective theory contains the same information as the corresponding one for QCD. However, while in 
QCD self-consistency is maintained for all momentum modes via the solution of the stationarity condition (|73|l . in 
the effective theory self-consistency is only required for the relevant quark and soft gluon modes. In this sense, the 
effective theory provides a simplification of the full problem. 



B. Dyson- Schwinger equations for relevant quarks and soft gluons 

After having specified r 2 in Eq. i|9U|) , we are now in the position to write down the Dyson-Schwinger equations 174|) 
explicitly. For the full inverse propagator of soft gluons we obtain with Eqs. I|74a|l . (|75a|l . (|85l) . and (|90|l 
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tv„ (g 0<22 rg 0<22 r +2Tr„ sr g 0<22 r +Tr„ (gvgr 



(91) 



The first term in square brackets takes into account the effect of quark-antiquark excitations as well as quark-hole 
excitations far from the Fermi surface. The second term is the contribution from excitations where one quark is close 
to the Fermi surface (a relevant quark) while the second is far from the Fermi surface or an antiquark (an irrelevant 
quark). The relevant quark propagator g can have diagonal elements in Nambu-Gor'kov space, corresponding to 
normal propagation of quasiparticles, as well as off-diagonal elements, corresponding to anomalous propagation of 
quasiparticles, cf. Eq. (|77|l . However, in the second term in square brackets the latter contribution is absent, because 
So, 22 is purely diagonal in Nambu-Gor'kov space, cf. Eq. H19[). This is different for the last term in square brackets, 
which corresponds to quark-hole excitations close to the Fermi surface. Both quark propagators have to be determined 
self-consistently and may have off-diagonal elements in Nambu-Gor'kov space. Consequently, the trace over Nambu- 
Gor'kov space gives two contributions, a loop where both quarks propagate normally, and another one where they 
propagate anomalously. Diagrams of this type have been evaluated in Ref. an( i l ea d to the Meissner effect for 
gluons in a color superconductor. 

For the full inverse propagator of relevant quarks we obtain with Eqs. (|74b|) . i|75b|) . l|89|) . and i|90|) 
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Tr 9 ArSo,2 2 r +Tr 9 A , 22 rgr +Tr 9 Argr 



(92) 



The first two terms in square brackets do not have off-diagonal components in Nambu-Gor'kov space. They contribute 
only to the regular quark self-energy. The other two terms in square brackets have both diagonal and off-diagonal 
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components in Nambu-Gor'kov space. The diagonal components contribute to the regular quark self-energy, in 
particular, the fourth term leads to the quark wave- function renormalization factor computed first in Ref. [33| . It 
gives rise to non- Fermi liquid behavior 11] . The off-diagonal components enter the gap equation for the color- 
superconducting gap parameter. 

The system of Eqs. (|91|) and 192|) has to be solved self-consistently for the full propagators of quarks and gluons. 
However, as was shown in Ref. |34| . in order to extract the color-superconducting gap parameter to subleading order it 
is sufficient to consider the gluon propagator in HDL approximation; corrections arising from the color-superconducting 
gap in the quasiparticle spectrum are of sub-subleading order in the gap equation. For our purpose this means that 
it is not necessary to self-consistently solve Eq. (|91|) together with Eq. (|92|l ; we may approximate Q on the right-hand 
side of Eq. (|91l) by Go,u- In essence, this is equivalent to considering only the first term on the right-hand side 
of Eq. I|92l) when solving Eq. I|91|l . Of course, under this approximation the effect of the regular quark self-energy 
(leading to wave-function renormalization) and of the anomalous quark self-energy (which accounts for the gap in the 
quasiparticle excitation spectrum) are neglected. 

With this approximation, and using Go = Go,ii © £0,22, we may combine the terms in Eq. Ij91|l to give 

-,2 



A aii + y Tr ?(^fg f) . (93) 



Taking the gluon cut-off scale A g to fulfill gfi <C A g < soft gluons are defined to have momenta of order gfi. We 
compute the fermion loop in Eq. (|93[) under this assumption (taking the soft gluon energy to be of the same order of 
magnitude as the gluon momentum). We then realize that the soft gluon propagator determined by Eq. (|93|l is just the 
gluon propagator in HDL approximation. We indicate this fact in the following by a subscript, A = Ajjdl- Armed with 
this (approximate) solution of the Dyson-Schwinger equation l|91[) we now proceed to solve Eq. I|92[) . We consider 
the two independent components [G + ] _1 + S + and <E> + in Nambu-Gor'kov space separately. Due to translational 
invariance, it is convenient to define [G+]- 1 (K,Q) = (l/T)[G+]- x (i^) S^ Q , Y,+ {K,Q) = (1/T)T,+(K) 8$ <Q , and 
using Eqs. {I7|), {HU, CHH, (HSU, w e obtain the Dyson-Schwinger equation for [G+] -1 + S + , 

[G+]- 1 (X) + S+(X) = [G+ U ]-\K) 

T 



ff 



- ]T { [A 0)22 C (K-Q) + [AhdlC (K-Q)} jfj,T a G Q , 22 (Q) l»T b 
Q 

£ { [Ao )22 C (K-Q) + [AhdlC (K-Q)} l,T a G + (Q) l„T b . (94) 
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V 

y 

Note that the first sum over Q runs over irrelevant quark momenta, < q < /j — A q and /1 + A q < q < 00, while the 
second sum runs over relevant quark momenta, fi — A q < q < /i + A q . There is no double counting of gluon exchange 
contributions, since the hard gluon propagator Ao,22 has support only for gluon momenta |k — q| > A g , while the 
HDL propagator is restricted to gluon momenta |k — q| < A g . To subleading order in the gap equation, we do not 
have to solve this Dyson-Schwinger equation self-consistently. It is sufficient to use the approximation G + — Gj 1:L on 
the right-hand side of Eq. (|94l) and to keep only the last term which, as discussed above, is responsible for non- Fermi 
liquid behavior in cold, dens e q uark matter. The net result is then simply a wave-function renormalization for the 
free quark propagator Gq U 33], 

M 2 

[G + ]-\K) + Z+(K) ~ [G+^-^K) + f k 70 In -rj- = [^ _1 (*o) fc + m] 70 - 7 ■ k , (95) 

fc 

where g = g/(3\/2V) and M 2 = (37r/4) m 2 , w ith the gluon mass parameter m g defined in Eq. (0. Neglecting effects 
from the finite life-time of quasi-particles 30] , which are of sub-subleading order in the gap equation, the wave- function 
renormalization factor is 

Z(ko)= (l+flnj^-j . (96) 

Due to translational invariance, it is convenient to define $ + (if, Q) = (1/T) & + (K) S^q and E + (K, Q) = 
TS + (K) 5^q, and the Dyson-Schwinger equation for <& + (K) becomes 

<f+(K) = .9 2 I E { t A ^2C {K-Q) + [AhdlC (K - Q) } l,(T a f 5+(Q) lu T» . (97) 
Q 
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Here, the sum runs only over relevant quark momenta, /i — A g < q < fi + A q . This is the gap equation for the 
color-superconducting gap parameter within our effective theory. There is no contribution from irrelevant fermions, 
since their propagator is diagonal in Nambu-Gor'kov space. 

While the gluon cut-off was taken to be A g < fj,, so that soft gluons have typical momenta of order gfi, so far 
we have not specified the magnitude of A q . In weak coupling, the color-superconducting gap function is strongly 
peaked around the Fermi surface 0,0, Q- For a subleading-order calculation of the gap parameter, it is therefore 
sufficient to consider as relevant quark modes those within a thin layer of width 2A q around the Fermi surface. For 
the following, our principal assumption is A q < g[i <C A g < //. As we shall see below, this assumption is crucial 
to identify sub-subleading corrections to the gap equation (j^IJ, which arise, for instance, from the pole of the gluon 
propagator. Note that this assumption is different from that of Refs. [llALjjj], where it is assumed that A q — A g . 

For a two-flavor color superconductor, the color-flavor-spin structure of the gap matrix is Q 

$+(K) = J 3 T 2l5 A+e(A q -\k-^\)<f>(K) , (98) 

where {J 3 )ij = — i&ijs and (T 2 )f g = —itfg represent the fact that quark pairs condense in the color- antitriplet, flavor- 
singlet channel. The Dirac matrix 75 restricts quark pairing to the even-parity channel (which is the preferred one 
due to the U(1)a anomaly of QCD). In the effective action 1|53JI . antiquark and irrelevant quark degrees of freedom 
are integrated out. The condensation of antiquark or irrelevant quark pairs, while in principle possible, is thus not 
taken into account; the bilocal source terms in Eq. (fTTf) only allow for the condensation of relevant quark degrees of 
freedom. The condensation of antiquarks or irrelevant quarks could also be accounted for, if one introduces bilocal 
source terms already in Eq. {Ejl , i.e., prior to integrating out any of the quark degrees of freedom. While there is 
in principle no obstacle in following this course of action, it is, however, not really necessary if one is interested in a 
calculation of the color-superconducting gap parameter to subleading order in weak coupling: antiquarks contribute 
to the gap equation beyond subleading order |3fij . and the gap function for quarks falls off rapidly away from the 
Fermi surface, i.e., in the region of irrelevant quark modes, and thus also contributes at most to sub-subleading order 
to the gap equation. Consequently, the Dirac structure of the gap matrix l|97[) contains only the projector A^ onto 
positive energy states. The theta function accounts for the fact that the gap function (j>(K ) pertains only to relevant 
quark modes. 

Inserting Eq. (|95|l and the corresponding one for [G~] _1 + X - , as well as Eq. (|98|l . into the definition l|79|l for the 
anomalous quark propagator, one obtains 



HQ) 

[qo/Z( qa W-et 



S+(Q) - J3T275 A q 0(A, - \q - ^ , J^/ 2 2 . (99) 



One now plugs this expression into the gap equation Ij97|l . multiplies both sides with J3T275Aj^, and traces over color, 
flavor, and Dirac degrees of freedom. These traces simplify considerably since both hard and HDL gluon propagators 
are diagonal in adjoint color space, [Ao,22]ah = $ab Aq22' [AuDiJah = ^afc A^p L . The result is an integral equation 
for the gap function </>(i"Q, 

K K ) = J y E [ A o^(K -Q) + A^ DL (K - Q) ] Tr s (A+7^7,) g2 . (100) 



Q 

The sum over Q runs only over relevant quark momenta, \q— fi\ < A q . Also, the 3-momentum k is relevant, \k— fi\ < A q . 

C. Solution of the gap equation 

In pure Coulomb gauge, both the hard gluon and the HDL propagators have the form 

A 00 (P) = A^(F) , A o,; ( J P) = , A ij (P) = (5 ij -p i f)A t (P) , (101) 
where A '* are the propagators for longitudinal and transverse gluon degrees of freedom. For hard gluons 

A e . 22 (P) - . ( 102a ) 
Al 22 (P) = , (102b) 
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while for soft, HDL-resummed gluons 
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(104a) 
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The HDL propagators l|103|l have quasiparticle poles at po = ±uoi^{p) 1 and a cut between po — —p and po = p [23l |. 
The gluon energy on the quasiparticle mass-shell is always larger than the gluon mass parameter, uji y t{p) > m g , where 
the equality holds for zero momentum, p — 0. 




(a) (b) 

FIG. 20: (a) The contour C in Eq. 11051 encloses the poles of tanh[go/(2T)] on the imaginary go axis, (b) Deforming the 
contour C and adding semicircles at infinity to enclose the poles of the quark propagator on the real qo axis. 

We hrst perform the Matsubara sum, using the method of contour integration in the complex go plane p3l l3lT | , 

t e /(«o ^^£ dq °l tanh (§0 /(go) (io5) 

where the contour C consists of circles running around the poles uj l n = (2n + 1)tvT of tanh[go/((2T)] on the imaginary 
qo axis, cf. Fig. [201(a). Inserting the propagators (|102J) and (|103[) into Eq. HlOOjl . we have to compute four distinct 
terms. The hrst one arises from the exchange of static electric hard gluons. Since Ag 2 2(P) does not depend on 
Po = ko — qo, only the quark propagator gives rise to a pole of f(qo), cf. Fig. [201 (b). After deforming the contour 
and closing it at infinity as shown in Fig. 1201 (b), one employs the residue theorem to pick up the poles of the quark 
propagator, 

TJ2 K.AP) F /7 f Q l 2 = i tanh (%) ^ [^,q) + ^(-S,,q)] , (106) 
V [qo/Z(qo)Y -e* p 2 \2T J 4e 9 

with e q = e q Z(e q ). Here, we have used the fact that the quark wave-function renormalization factor is an even function 
of its argument, Z(qo) = Z(—qo), cf. Eq. JOEJl. An essential assumption in order to derive Eq. p()6|l is that the gap 
function <p(Q) is analytic in the complex go plane. This assumption will also be made in all subsequent considerations. 

By the same method one computes the second term in Eq. (|100|) . corresponding to magnetic hard gluon exchange. 
This is slightly more complicated, since not only the quark propagator but also A 22 (P) has poles at p = ±p, which 
are located at qo = k ± p in the complex q plane, cf. Fig. 1211 The external quark energy k Q is hxed and, prior to 
analytic continuation fcp ~> £fc + *?7 to the quasiparticle mass-shell, is equal to one particular fermionic Matsubara 



27 




FIG. 21: Same as in Fig. 1201 but for magnetic hard gluon exchange. Now also the gluon propagator has poles at ko ± p in the 
complex go plane. These are further away from the imaginary axis than the poles e q of the quark propagator, because for our 
choice of quark and gluon cut-offs, A, «A 9 , we have e, < A, C A 9 < p. 



frequency, cf. Fig. 1211 The residue theorem now yields four contributions, two from the quark and two from the gluon 
poles. Using tanh[(fco ±p)/(2T)] = ±coth(p/2T) and analytically continuing fcp — > e k + if] we find 
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.(107) 



Since the gluon momentum is hard, p > A g , and thus much larger than the quasiparticle energies e k , i q which are at 



most of the order of the quark cut-off A g <C A s , to order 0(A q /A g ) we may neglect the terms (e k ±e q 
denominators of the first term. Furthermore, in the second term we may approximate Z(p ± ik) 



-irf) 2 in the energy 
Z(p) = 1 + 0(g 2 ) 



and 4>{p ± £fc,q) — <j>(p,q). Note that the gap function is far off-shell for p > A g ^> A q > \fj, — q\. Then, to order 
0(A q /Ag), we may also neglect e k , e q in the energy denominators of the second term. We obtain 
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Let us estimate to which order the two remaining terms contribute to the gap equation 1100|) . At T = 0, we may set 
the hyperbolic functions to one. We shall also ignore the difference between the on-shell and off-shell gap functions, 



and take (f>(p, q) 



(±e 9 ,q) 



= const.. For the purpose of power counting, we may restrict ourselves to the 



leading contribution of the Dirac traces in Eq. (|100f> . which is of order one, cf. Eqs. 
the leading contribution of the first term in Eq. I|l()8|l . we may also set Z 2 (e q ) ~ 1 
magnitude of the quark momentum is J dqq 2 , while the angular integration is J dcos9 
term in Eq. (|108|l leads to the following contribution in the gap equation 
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(109) 



where we approximated k 



a and employed the weak-coupling solution J5J to estimate ln(2A g /0) ~ 1/g. 



Furthermore, for A g < //, the angular logarithm is ln(2/x/A s ) ~ 0(1). According to the discussion presented in the 
introduction, the contribution from hard magnetic gluon exchange is thus of subleading order in the gap equation. 
Note that the term arising from hard electric gluon exchange, Eq. I|106|) . is °f the same order as the first term in Eq. 
(|108fl . and thus also contributes to subleading order. The way we estimated the first term on the right-hand side of 
Eq. (|108J) is equivalent to just taking the hard magnetic gluon propagator in the static limit, Aq 2 2(P) — 1/p 2 j which 
is correct up to terms of order 0(A 2 /A 2 ). To this order, the propagator for hard magnetic gluons is thus (up to a 
sign) identical to the one for hard electric gluons. Since the ratio A q /A g ~ g/i/ /i = g, this approximation introduces 
corrections at order 0(g 3 4>) in the gap equation, which is beyond sub-subleading order, 0(g 2 (j>). 
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Similarly, we estimate the contribution of the second term in Eq. I|108[) to the gap equation (|100fl . 



9 



dp 2 A q 3 
dqq —~g4>- g q 

A a P 2 A 9 



(110) 



i.e., for our choice A q /A g ~ g, this term contributes beyond sub-subleading order. Note that this estimate is conser- 
vative, as we assumed the off-shell gap function to be of the same order as the gap at the Fermi surface, <fr(p, q) ~ <j>. 
However, we know Q that, for energies far from the Fermi surface, e q ~ A q < gfi : even the on-shell gap function is 
suppressed by one power of g compared to the value of the gap at the Fermi surface, 0(A g ,q) ~ g<p. The off-shell 
gap function at go = P ^ A g ^> A q may be even smaller. In order to decide this issue, one would have to perform a 
computation of the gap function for arbitrary values of the ener gy q n , and not just on the quasiparticle mass-shell, 
go = e q . We note that for the choice A g ~ A g for the cut-offs fla Il7|. the ratio A q /A g is of order one and cannot be 
used as a parameter to sort the various contributions according to their order of magnitude. The expansion of the 
denominators in powers of A q /A g as seen on the right-hand side of Eq. H108|) is then inapplicable. 
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FIG. 22: Evaluating the Matsubara sum for HDL-resummed gluon propagators, (a) The original contour C in Eq. 11051 . There 
is no circle around the point fco = go, where the corresponding term in the Matsubara sum has a cut arising from the HDL 
gluon propagator, (b) Deforming the contour C. (c) The contour C cut running around the cut. (d) The contour C\ = C + C cu t 
which is closed at infinity. 

The third and fourth terms in the gap equation H1UU|) arise from soft, HDL-resummcd electric and magnetic gluon 
exchange. Evaluating the Matsubara sum via contour integration in the complex go plane is considerably more difficult 
than in the previous cases, because the HDL gluon propagators A^jp L do not only have poles but also cuts. The 
analytic structure is shown in Fig. 1221(a)- Besides the poles of the quark propagator at go = ±e g , there are also those 
from the gluon propagator at go = fco ±W£.t(p). The cut of the gluon propagator between — p < pg < p translates into 
a cut between ko — p < qo < ko + p. Prior to analytic continuation, the gluon poles and the cut are shifted away from 
the real axis and located at the (imaginary) external Matsubara frequency fco- 

The Matsubara sum over go is evaluated in the standard way, cf. Eq. (|105f) . with the caveat that the contribution at 
go = fco, where the cut of the gluon propagator is located, has to be omitted. This is similar to the zero-temperature 
case where the Matsubara sum becomes a continuous integral along the imaginary go axis and where one has to avoid 
integrating over the cut. Alternatively, the term go = fco can be included in the Matsubara sum if one shifts the cut by 
some small amount ±ie along the imaginary go axis. The final result will be the same, as one still has to circumvent 
the cut by a proper choice of the integration contour. 
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We now deform the contour as shown in Fig. [221 (b), and add and subtract a contour integral running around the 
cut, Fig. [221 (c). The integral over the contour C + C cut can be closed at infinity, yielding the contour C\ shown in Fig. 
[22(d). One obtains 
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Evaluating the integral over C\ is rather similar to the case of hard gluon exchange: one just picks up the poles of the 
quark and gluon propagators inside the contour C\ . After analytic continuation fcg — ► £fc + if] one obtains 
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Here, we approximated the quark wave-function renormalization factor Z(u>i t t ± £&) ~ 1 + 0(g 2 ). We also expanded 

[1 + 0(e^/w| f )]. For our choice of the cut-off 



the denominators of the quark propagator (?& ± u>t t t + iff) ~ e q — UJ it\ 
Aq ^ gfJ. ~ m 9 , we may estimate > m 9 > A g > e q , i.e., the corrections of order 0(e 2 /u> 2 t ) are small everywhere 
except for a small region of phase space where p ~ and e q ~ A 9 . (In principle, in the expansion of the denominators 
there are also linear terms, ~ ±ik/toe,t, but these are very small everywhere for external momenta close to the Fermi 
surface, k ~ fi.) Note that the gap function is again off-shell at the gluon pole, although not as far as in the case of 
hard gluon exchange, cf. Eq. (|108fl . The residues of the HDL gluon propagators at the respective poles are [23| 
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To very good approximation, one finds that Z t (oj t ,p) — — l/(2w t ) for all momenta p. In the longitudinal case, the 
residue is very well approximated by Zi{u>i,p) ~ —u>i/(2p 2 ) for small momenta p < m g , while for large momenta, 



<p, Zi(u>i,p) ~ exp[— 2p 2 /(3m 2 )]/p 1 i.e., it is exponentially suppressed |36j . 



These approximate forms allow for a simple power counting of the gluon-pole contribution in Eq. (|112fl to the gap 
equation (|100|) . To this end, we approximate the gap function by its value at the Fermi surface, 0(±w^ + efc, q) — 4>, 
and consider the limiting case T = where coth[wf ! t/(2T)] = 1. Then, the contribution from the longitudinal gluon 
pole is 
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In the first p integral, which only runs up to the scale m s , one may approximate lui ~ m gi while in the second p 
integral, which runs from m g to A g < /i, one may take uj£ ~ p. To obtain the right-hand side of Eq. I|114|) we have 
set k ~ q ~ p,, and we have employed our choice A q < g[i for the quark cut-off. This also allowed us to approximate 
logarithms of A q /m g by numbers of order 0(1). With this choice for the quark cut-off, the contribution l|114|) is of 
sub-subleading order, ~ 0(g 2 (j)), to the gap equation. 

With a more careful evaluation of the integrals, one could extract the precise numerical prefactor of the sub- 
subleading contribution (f 1 1 4|) . Note, however, that further suppression factors may arise from the off-shcllness of the 
gap function at ifi^toej + q), which consequently would render this contribution beyond sub-subleading order. As 
noted previously, this issue can only be decided if 4>(qo, q) is known also off the quasiparticle mass-shell, and not only 
on-shell. We also note that the 1/p 2 factor in the residue Zg is an artifact of the Coulomb gauge |36(, and does not 
appear in e.g. covariant gauge. One would have to collect all other terms of sub-subleading order to make sure that 
the complete sub-subleading contribution is gauge invariant and the term (|114f> not cancelled by some other terms. 

Similarly, we estimate the contribution from the transverse gluon pole, 
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where we defined £ = q — fx. We approximated dpp ~ dujt u)t since, for the purpose of power counting, to very good 
approximation one may take the dispersion relation of the transverse gluon equal to that of a relativistic particle with 
mass m g , uj t (p) ~ (p 2 + m 2 ,) 1 / 2 . We also used A q < m g <C A g and fc ~ q ~ /i. In conclusion, also the transverse 
gluon pole possibly contributes to sub-subleading order in the gap equation, with the same caveats concerning the 
off-shellness of the gap function as mentioned previously. 

Let us now focus on the integral around the cut of the gluon propagator in Eq. If) . We substitute qo by 
Po = fco — qa = uj and use the fact that tanh[qo/(2T)] = — coth[w/(2T)]. Since the gluon propagator is the only part 
of the integrand which is discontinuous across the cut, we obtain after analytic continuation fco — > £fc + if] 

-2S L \ tanh (I) A ™ l(P) [qj{q%-,\ 

fP , 1 , , f u \ Z 2 (i k -uj)(t){e k -uj,q) e t 

dw-coth — — p—r r — ^ p'Juj, p) , (116) 

, 2 \2TJ (e k -uj + ir]) 2 -{Z(e k -uj)e q } 2 cut 

where p^ut {oj,p) = ImA^jp L (uj + irj, p) /n is the spectral density of the HDL propagator arising from the cut. Explicitly, 
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In order to power count the contribution from the cut of Afj DL to the gap equation, it is sufficient to approximate 
the spectral density by Q 
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This form reproduces the correct behavior for uj <C p. For uj < p, it overestimates the spectral density when p 



< 



rn„ 



while it slightly underestimates it for p 



> 



For the gap equation, however, this region is unimportant, since the 



respective contribution is suppressed by the large energy denominator (e k — uj + irj) 2 — [Z(e k — uj) e q ] 2 ~ p 2 in Eq. 
(|116fl . To leading order, we may set Z(e k — uj) ~ 1. We also approximate 4>(e k — uj, q) ~ cf>. Then, the uj integral can be 
performed analytically. (One may compute this integral with the principal value prescription; the contribution from 
the complex pole contributes to the imaginary part of the gap function, which we neglect throughout this computation.) 
This produces at most logarithmic singularities, which are integrable. We therefore simply approximate the uj integral 
by a number of order 0(1). Consequently, the contribution from Eq. (|116fl to the gap equation is of order 




(119) 



< 



where we approximated the p integral by a method similar to the one employed in Eq. 1)1140 . For our choice A q 
Eq. i|119D constitutes another (potential) contribution of sub-subleading order to the gap equation. 

Finally, we estimate the contribution from the cut of the transverse gluon propagator. For all momenta p and 
energies — p < uj < p, a very good approximation for the spectral density (|117bfl is given by the formula 



Pcut(^:P) 



M 2 



uip 



P 6 + (M 2 uj) 2 



(120) 



This approximate result constitutes an upper bound for the full result (|117bjl • The advantage of using this approximate 
form is that, interchanging the order of the p and uj integration in the gap equation, the former may immediately be 
performed. Approximating Z{e k — uj) ~ 1, neglecting the dependence of the gap function on the direction of q, and 
defining A = max(|fc — q\,uj), at T = the contribution to the gap equation is 



(121) 



1 




e k + uj + e q 
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Here, we have used the fact that the particular combination of arctan's in the first line effectively cuts off the u> 
integral at the scale u> ~ M. As usual, we have set k ~ q ~ fi. If we simply neglect the off-shell behavior of the 
gap function and approximate 4>(ik ± w, <?) — <j>: this contribution would (at least) be of subleading order. Note that 
the corresponding contribution in previous treatments of the QCD gap equation, cf. for instance Eq. (67) of Ref. ||, 
was discarded as being of higher order. At this point, we refrain from a more careful evaluation of the contribution 
1)121(1 . because this requires a calculation of the gap function off the quasiparticle mass-shell. Since the purpose of 
the present work is to show that our method reproduces previous results, we follow Ref. [j| and also discard the 
contribution l|121|) in the following. 

The remaining term from the evaluation of the Matsubara sum in Eq. I|lll|) is the contribution from the quark 
pole, i.e., the first line of Eq. 1)1120 . This has to be combined with the subleading-order terms from hard-gluon 
exchange, i.e., from Eq. l|lUtj[) and from the first line of Eq. (|108fl . in order to obtain the gap equation which contains 
all contributions of leading and subleading order. Before doing so, however, we also evaluate the Dirac traces in Eq. 
(|100[> . In pure Coulomb gauge, we only require 



(k + q) 2 - p 2 
2kq 

p 2 (k 2 - q 2 ^ 2 



Tr s (A+ 7o A q7o ) = P , (122a) 



^ - M = 2 - fk-q + -Ykq^T ■ (122b ) 

where we used p 2 = (k — q) 2 = k 2 + q 2 — 2 k gk • q to eliminate k • q in favor of p 2 . Let us estimate the order of 
magnitude of the terms arising from the traces at the Fermi surface, k = fx. Setting q = fi + £, where — A q < £ < A 9 , 
one obtains 

Tr s (A+ 70 A q 70 ) = 2 - ^- + O {^j , (123a) 

2 / C2 ■ 



(<F - m Tr s (A+ 7i A- 7j ) = -2 - ±- + O (L^j . (123b) 



As shown above, the contribution from hard-gluon exchange is at most of subleading order. Thus, for this contribution 
it is sufficient to keep only the leading terms in Eq. (|123|l . i.e., one may safely neglect terms of order 0(£, 2 /fi 2 ) < 
0(A 2 /A 2 ) ~ 0(g 2 ) or higher. Note that, since for hard gluon exchange p ~ /i > A ff , the terms p 2 /(2kq) cannot be 
omitted. However, since the magnetic gluon propagator is effectively ~ 1/p 2 , cf. Eq. (|108|l . i.e., (up to a sign) identical 
to the electric propagator, these terms will ultimately cancel between the electric and the magnetic contribution. This 
cancellation is well-known, see for instance Ref. |33, and is special to the spin-zero case. It does not occur in spin- 
one color superconductors where there is an additional exponential prefactor which suppresses the magnitude of the 
spin-one gap relative to the spin-zero case [37) . 

As is well-known, electric soft-gluon exchange also contributes to subleading order in the gap equation. Thus, as 
in the case of hard-gluon exchange, we may drop the terms of order 0(£ 2 //i 2 ) in Eq. I|123a|l . On the other hand, 
magnetic soft-gluon exchange constitutes the leading order contribution to the gap equation. We therefore would 
have to keep all terms up to sit&leading order, i.e., ~ 0(£//i). Fortunately, the corrections to the result (|123b)l are of 
order 0(£ 2 //x 2 ) ~ 0(g 2 ), i.e., they are of sit&-sit&leading order and thus can also be omitted. 

We combine Eqs. <|106[1 . (|108|) . and the first line of Eq. (|112|l . and assume that the gap function is even in its energy 
argument, <j>{— e g ,q) = (p(e q ,q), and isotropic, <p(e q ,q) = <fi(e q ,q) = <fi q . Then, on the quasiparticle mass-shell fco = ?fc 
the gap equation I1UUI) becomes 
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(124) 



The next step is to divide the integration region in the p — q plane into two parts, separated by the gluon "light cone" 
| efc — se q \ = p. For our choice A q <C A g the region, where \ik — se q \ < p, is very large, while its complement is rather 
small. In order to estimate the contribution from the latter to the gap equation, we may approximate the HDL gluon 
propagators by their limiting forms for large gluon energies, cf. Eqs. (|103fl . (|104fl . 



Po > P 



(125) 



.32 



Following the power-counting scheme employed previously, the contribution from the electric sector is of order 

c\lk-se q \ 

P 



a 4> f +A * q (e fc - 8 g g ) 2 

K Jfj.-Ag e q m g J\k-q\ 



dp 2 <i> f A " , f 2 , A 2 



"g JQ 



9 



(126) 



This is a contribution of sub-subleading order, as long as one adheres to the choice A 9 < g/i. Analogously, we estimate 
the contribution from the magnetic sector to be 
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dq 
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Consequently, all contributions from the region — se q \ > p are of sub-subleading order, and the further analysis 
can be restricted to the region \ik — se q \ < p. In this region, it is permissible to use the low-energy limit of the HDL 
gluon propagator, which follows from Eqs. (|103fl . <|1U4[1 keeping only the leading terms in the gluon energy, 



Po < P ■ 
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(128) 



Here, we only retained the real part of the transverse gluon propagator, since the imaginary part contributes to the 
imaginary part of the gap function, which is usually ignored. (In Ref. [j| it was argued that, at least close to the 
Fermi surface, the contribution of the imaginary part is of sub-subleading order in the gap equation.) With the 
approximation 1)128(1 . the gap equation (|124fl becomes 
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where we already performed the integration over hard gluon momenta p > A g . The integration over soft gluon mo- 
menta can also be performed analytically. Formally, the terms ~ p 2 /(4kq) give rise to subleading-order contributions, 
~ A 2 /(8fcg), but they ultimately cancel, since they come with different signs in the electric and the magnetic part. 
Other contributions from these terms are at most of sub-subleading order. Exploiting the symmetry of the integrand 
around the Fermi surface and setting k ~ /i, we arrive at 
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(130) 



Here, we have neglected terms ~ £k — se q against 3 m 2 under the logarithm arising from soft electric gluons, and terms 
~ (ik — se q ) 6 against M 4 (ik — se q ) 2 under the logarithm from soft magnetic gluons. 



Now observe that the gluon cut-off A g cancels in the final result, 



18tt 2 



d(q - fi) 



Z\e q ) 



tanh 
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where b = 256n 4 [2/(N f g 2 )} 5 ^ . This is Eq. (19) of Ref. [T^, 
integration, 8, replaced by the quark cut-off A q . 



since g 2 = g 2 /(187r 2 ), with the upper limit of the (q — fx) 



The solution of the gap equation p31|l is well-known, and given by Eq. J5J. As was shown in Ref. @, the dependence 
on the cut-off A q enters only at sub-subleading order, i.e., it constitutes an O(g) correction to the prefactor in Eq. 
(J2J). Therefore, to subleading order we do not need a matching calculation to eliminate A q . 

The result (f 1 311) shows that the standard gap equation of QCD can be obtained from the effective action 115311 . The 
above, rather elaborate derivation of Eq. (|131|) demonstrates that, in order to obtain this result, it is mandatory to 
choose A ? < A 9 . This also enabled us to identify potential sub-subleading order contributions. However, we argued 
that, at this order, the off-shell behavior of the gap function has to be taken into account. 
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V. SUMMARY AND OUTLOOK 

In this paper we have presented a formal derivation of an effective action for non-Abelian gauge theories, Eq. I|53[) . 
We first introduced cut-offs in momentum space for quarks, A q , and gluons, A g . These cut-offs separate relevant from 
irrelevant quark modes and soft from hard gluon modes. We then explicitly integrated out irrelevant quark and hard 
gluon modes. The effective action (|53|) is completely general and, as shown explicitly in Sec. lIII A1 after appropriately 
choosing A q and A s , it comprises well-known effective actions as special cases, for instance, the "Hard Thermal Loop" 
(HTL) and "Hard Dense Loop" (HDL) effective actions. We also demonstrated, cf. Sec. IIII Bl that the high-density 
effective theory introduced by Hong and others 0, 0, 0, 0, 0, is contained in the effective action ifCT) . 

We then showed how the QCD gap equation can be derived from the effective action (15311 . The gap equation is a 
Dyson-Schwinger equation for the anomalous part of the quark self-energy. It has to be solved self-consistently, which 
is feasible only after truncating the set of all possible diagrams contributing to the Dyson-Schwinger equation. Such 
truncations can be derived in a systematic way within the Cornwall- Jackiw-Tomboulis (CJT) formalism p^j . Here, 
we only include diagrams of the sunset-type, cf. Fig. in the CJT effective action, which gives rise to one-loop 
diagrams (with self-consistently determined quark and gluon propagators) in the quark and gluon self-energies. 

Usually, the advantage of an effective theory is that the degree of importance of various operators can be estimated 
(via power counting) at the level of the effective action, i.e., prior to the actual calculation of a physical quantity. This 
tremendously simplifies the computation of quantities which are accessible within a perturbative framework. On the 
other hand, the requirement of self-consistency for the solution of the Dyson-Schwinger equation invalidates any such 
power-counting scheme on the level of the effective action. For instance, perturbatively, the right-hand side of the gap 
equation Q is proportional to g 2 . However, self-consistency generates additional large logarithms ~ ln(///0) ~ 1/g 
which cancel powers of g. 

Nevertheless, there is still a distinct advantage in using an effective action for the derivation and the solu- 
tion of Dyson-Schwinger equations for quantities which have to be determined self-consistently, such as the color- 
superconducting gap function in QCD. This advantage originates from the introduction of the cut-offs which separate 
various regions in momentum space. They allow for a rigorous power counting of different contributions to the 
Dyson-Schwinger equation. We explicitly demonstrated this in Sec. IIVI where we reviewed the calculation of the 
color-superconducting gap parameter to subleading order. 

In order to obtain the standard result (GJ), it was mandatory to choose A q < g/i <C A g < [i. This is in contrast 
to previous statements in the literature |l5l ll& Il7| that a consistent power-counting scheme requires A q ~ A g . In 
particular, the choice A q -C A g has the consequence that the gluon energy in the QCD gap equation is restricted to 
values po < A q , while the gluon momentum can be much larger, p < A g . This naturally explains why it is permissible 
to use the low-energy limit (II28[) of the HDL gluon propagators in order to extract the dominant contribution to the 
gap equation (which arises from soft magnetic gluons). In previous calculations of the gap within the framework of 
an effective theory 0,E3,E1> the low-energy limit for the HDL propagators was used without further justification, 
even though for the choice A g ~ A g the gluon energy can be of the same order as the gluon momentum. The physical 
picture which arises from the choice A q < g/i -C A g < ji is summarized in Fig. If 71 Relevant quarks are located within 
a thin layer of width ~ A q around the Fermi surface. Soft gluon exchange mediates between quarks within a "patch" 
of size ~ A g inside this layer. The area of the patch is much larger than its thickness. Hard gluon exchange mediates 
between quark states inside and outside of the patch. 

In the course of the calculation, we were able to identify various potential contributions of sub-subleading order. 
However, we argued that, at this order, a solution of the gap equation must take into account the off-shell behavior 
of the gap function. For a complete sub-subleading order calculation it also appears to be necessary to include 2PI 
diagrams beyond those of sunset topology in T2 , cf. Eq. (|90|l and Fig. ^| Besides an improvement of the result for 
the color-superconducting gap parameter beyond subleading order, we believe that our rather general effective action 
(15311 can serve as a convenient starting point to investigate other interesting problems pertaining to hot and/or dense 
quark matter. 
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